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SECTION  SEVEN-NUMERICAL 
.  METHODS  OF  SOLUTION 

USE  OF  NUMERICAL  ANALYSIS  IN  INTERNAL  BALLISTICS 

^Various  variable  quantities  possessing  definite  physical 

V“**rs'*  **'--'*'  , 

significance  usually  take  part  in  processes  which  occur  in  nature 
or  are  conside red  in  technology.  In  this  connection,  numerical 
variations  of  one  or  more  quantities  are  accompanied  by  or 
.  .associated  with  variations  of  other  variables.  Thus  there  always 

'.exists  a  definite  functional  relation  between  the  variable  quantities 

V  . 

vnder  consideration.  This  functional  relation  may  be  expressed 
by  means  of  three  methods,  such  as  tabulations,  diagrams,  and 

^  formulas.  In  .the  vast  majority  of  processes,  especially  those 

* 

encountered  in  technology,  this  relation  is  expressed  by  the  aid 

•  'f 

of  {tables  or  diagrams  obtained  directly  from  experiment  or  from 
observation  of  the  process,  whereas  the  formulas  appear  only  after 
subsequent  analysis  of  the  results  obtained,  and  then  only  in  the 
case  of  the  simplest  processes.  It  is  thus  apparent  that  the  most 
natural  means  of  expressing  a  functional  relation  between  variable 
quantities  representing  in  their  totality  the  process  under 

investigation  is  a  tabulation;  this  is  especially  true  when  such 

r 

•’  a  process  is  being  utilized  directly  for  technological  purposes, 

Z_ln  which  case  a  formula  or  even  a  diagram  will  not  serve  the  purpose, 
:and  only  numerical  values  of  the  variable  quantity  considered  to 
be  of  primary  importance  on  the  basis  of  practical  considerations 
-  must  be  had. 

- Numerical  analysis  must  thus  serve  as  a  means  for  studying  and 

making  practical  utilization  of  the  functional  relations  between  the 
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variable  quantities  involved. 

Ordinary  differential  equations  may  bo  approximately  integrated 
by  means  of  any  of  the  known  methods  of  which  there  are  a  great  many. 
These  methods  include  the  following:  ‘  3— 

1)  Expansion  in  a  Taylor's  series  in  powers  of  the  argument. 

2)  Integration  by  the  method  of  successive  approximations. 

3)  Expansion  in  a  series  in  powers  of  small  parameters  entering 

into  the  equation.  _ 

4)  Expansion  in  a  series  in  powers  of  the  initial  values  of  the 
unknown  function  and  its  derivatives. 

5)  Method  of  successive  approximations  applied  to  equations  for 

1 

vibrational  motion. 

6)  Methods  of  Euler.,  E.L.  Bravin,  and  others. 

7)  Method  of  numerical  integration. 

The  first  six  methods  do  not  require  the  use  of  finite  differences; 
all  variants  of  the  method  of  numerical  integration  are  based  on  the 
use  of  such  differences. 

The  principal  variants  of  the  method  of  numerical  integration 
are  discussed  in  the  book  by  Academician  A.N.  Krylov,  "Lectures  on 
Approximate  Computations",  /~4_7 

Fundamental  information  on  the  theory  of  finite  differences  and 
on  the  technical  features  of  their  application  to  the  engineering  of 
artillery  may  be  found  in  the  book  by  Professor  G.V.  Oppokov,  "Numerical 
Analysis  Applied  to  the  Science  of  Artillery".  /"*5_7 

Internal  ballistics  is  an  applied  science  possessing  a 
perfectly  definite  technical  content  (the  study  of  the  motion  of 
a  projectile  in  the  bore  of  a  gun  and  of  the  laws  of  burning  of  powder) 
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and  a  perfectly  well-defined  technical  objective:  the  cx-eation  of 
eeans  for  plotting  the  curve  of  the  speed  of  the  projectile  in'' the 
bore. and  for  plotting  the  curve  of  the  pressure  of  the  powder 
gases  in  "the  bore  as  functions  of  the  path  of  the  projectile  and 


of  time. 


i  t^b' 


These  curves  can  be  plotted  after  obtaining  suitable  tabulations 
giving  the  functional  relation  between  the  various  variable 
quantities  participating  in  the  phenomenon  of  a  shot.  The  necessary 
tabulations  are  obtained  by  analyzing  primary  experimental  data  and 
those  formulas  which,  in  their  simplest  form,  express  the  relation 
existing  between  the  initial  variable  quantities. 

The  use  of  numerical  analysis  constitutes  the  subject  treate'd 

in  this  section  of  the  book. .  The  essence  of  numerical  analysis, 

its  specific  features  and  its  principal  operations  will  also  be 

considered  here  in  the  proper  degree.  _ _ 

CHAPTER!  -  NUMERICAL  INTEGRATION  BY  FINITE  DIFFERENCES 

(Written  .by  Professor  G.V.  Oppokov) 

1.  APPLICATION  OF  NUMERICAL  INTEGRATION  TO  THE ’ DETERMINATION 

OF  FUNCTIONS  ------ 


1)  Concept  of  Tabular  Functions  ;■ 

That  variable  function  whose,  numerical  values  can  be  chosen 
arbitrarily  is  usually  designated  as  the  argument  or  as  the 
independent  variable  quantity.  The  remaining  variable  quantities 
taking  part "in  the  process  under  consideration  are  then  designated 

as  funciTofiir  ...”  .  .......  ~ 

Let  us  assume  that  a  certain  independent  variable  quantity  takes 
on  a  series  of  particular  values: 
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which  arc  separated  from  one  another  by  invariable  equal  intervals,, 
so  that: 


*1  “  X0  "  X2  -  X1  "•••"  Xi+1  "  Xi  “•••"  Xn  -  Xn-1 


-  h. 


This  interval  is  called  the  step  of  the  argument  and  is 
designated  by  h. 

In  addition  to  the  step  h,  the  limits  of  the  region  x^  and  xn 
must  be  stated  in  the  form  of  finite  numerical  values. 

Cases  in  which  the  step  h  is  variable,  or  in  which  the  variable 
x  assumes  infinitely  large  values  in  the  region  under  consideration 
from  Xq  to  xn,  including  the  limits  of  the  region  themselves,  are 
not  considered  at  all. 

Furthermore,  let  some  other  variable  quantity  y  also  assume 
a  series  of  particular  values: 

V  yr  V'*”  yi'-*-'  v 

each  of  these  particular  values  corresponding  to  one  of  the  particular 
values  of  the  argument  x,  so  that  it  is  always  true_lhat: 


yi  "  f<xi>'  ‘ 


where: 

i  —  Op  1,  2, ,  n  1,  n. 
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It  follows  that-  the  functional  relation  between  the  variables 
y- and  x  is  established  in  the  form  of  a  tabulation. 

*  r  4  -  * 

An  example  of  such  a  relation  is  presented  in  Table  4. 

*•*  “  **’  *  / 

table  4  -  Pressure.  Curve  as  a  Function  of  Time  t. 


23 

23.5 

— 

24 

24.5 

25 

175 

223 

274 

330 

394 

0 

30.5 

30.75 

31 

1743 

1983 

2097 

2175 

- 

26  26.527.027.5 

546  636  ^39  857 


.  This  table  is  the  result  of  analysis  of  experimental  data 
obtained  by  burning  in  a  bomb  a  weighed  sample  of  strip  powder  grade 
SP  (1  x!8  x  40  am)  at  A-  0.201.  The  values  of  t  were  chosen  when 
evaluating,  the  experimental  data,  so  that  t  is  the  argument;  the 
values  of  pressure  p  were  taken  from  the  experimental  curve  for  the 
chosen  values  of  t.  Consequently,  the  pressure  of  the  powder  gases 
is  a  function  of  the  time  t. 

-“  Use  has  been  made  in  Table  4  of  the  so-called  vertical  notation 
which  we  shall  adopt  henceforth.  In  this  notation  the  particular 
values  of  each  of  the  variable  quantities  are  invariably  placed  in 
one  row,'  while  each  column  contains  one  particular  value  of  each  of 
these  variables.  This  notation  is  found  to  be  most  convenient  for 
the  various  operations  to  which  the  functions  stated  by  the 
tabulation  are  subjected. 

2)  Finite  Differences  of  Various  Orders. 


Once  a  function  has  been  given  in  table  form,  it  is  not  difficult 

,  '  .  / 

~  ~to  find  the  so-called  finite  differences,  of  this  function  or,  more 
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simply,  its  table  of  differences,  these  differences  being  of  various 
orders.  Thus,  for  a  portion  of  Table  4,  we  obtain  the  following 
table  of  differences. 


Table  4-a. 


t-103 

23 

23.5 

'24 

24.5 

25 

25.5 

p 

175 

223  • 

274 

330 

:i  ■  i 

466 

Ap 

48 

51 

56 

64 

72’ 

- 

A3p 

3 

5 

8 

8 

- 

- 

A3p 

2 

3 

0 

- 

- 

- 

Thus  on  the  basis  of  the  table  of  particular  function  values 
it  is  possible  to  compute  the  following  differences: 

Ap0  «  Px.  -  PQ  ~  223  -  175  -  48 
Ap^  *>  Pg  ~  P^  "  274  -  223’  «  51 
Ap2  -  P3  -  P2  -  330  -  274  -  56 


(?•>' 


o 


These  new.  differences  are  now  called  differences  of  the  second 
order,  or  second  differences. 

•  From  the  second  differences  it  is  also  possible  to  compute  in 
a  similar  manner  the  third  differences: 

A3P0  -  A2Pi  -  A2P0  -  5  -  3  -  2 
A3?^  -  A2p^  -  A2Pj  -  8  -  5  -  3 
A3P2  “  A3p<j  -  A2P2  **8-8-0 


followed  by  the  fourth  differences,  etc.  Generally,  by  definition, 
the  kth  difference  equals: 

Ai»l  •aW’i‘ 

6  Rule.  In  formulating  a  table  of  differences,  the  number  at 

the  leftmost  be  subtracted  from  the  number  at  the  right  in  the 

«- 

same  row,  and  the  result  recorded  in  the  next  lower  row  under  the 
.  number  at  the  left. 

It  is  useful  to  point  out  that  in  following  this  rule  the 
differences  with  the  same  subscript  (the  function  itself  also  being 
considered  as  a  zero  order  difference)  are  automatically  recorded 
in  one  and  the  same  column. 

•  "  ,  W 
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Table  4-b. 


X 

*o 

X1 

X2 

X3 

X4 

y 

yo 

y i 

*2 

y3 

y4 

Ay 

Ayo 

Ayi 

Ay2 

Ay3. 

- 

A2y 

&2y0 

&2/i 

A2y2 

- 

A3y 

A3y0 

A3y! 

*** 

“ * 

In  the  table  of  differences,  in  the  column  for  i  -  0,  for 
example,  are  found  differences  such  as: 


48, 

2 


etc. 

3)  Certain  Properties  of  Finite  Differences. 

1.  All  differences  of  any  order  can  be  expressed  only  in  terms 
■  of  particular  functions  of  the  .function  itself. 

By  definition  of  the  first  differences: 

A*0  '  yl  '  y0'  Ayl  "  y2  -  *1’  ■  y3  "  y2* 

Then,  by  definition  of  the  second  differences: 

A\  -Ayj  -Ay0  -  <y2  -  yj)  -  (yx  -  y0>  -  y2  -  2y,  +  y0: 

* 

A2yx  -Ay2  -  Ayx  -  (y3  -  y2>-  (y2  -  y^  -  y3  -  2y2  +  yx 

etc.  Generally: 
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A  yi  -  yi+2  -  ?yi+i  +  >v 


The  same  procedure  can  also  be  applied  to  differences  of  higher 


orders.  ■  _ - 

2.  A  constant  number  can  be  taken  outside  the  difference  symbol 
of  any  order  and,  conversely,  can  be  brought  inside  this  difference 
symbol.  >' 

Let : 


y  -  c  •  f (x) . 


where  c  is  a  constant  number.  On  the  basis  of  the  definition  of 
differences, -  we  have: 


Ay  -  c*f (x  +  h)  -  c*  f (x)  *  c-/~ f (x  +  h)  -  f (x)  7  -  c-  A f (x) . 


Thus, 


A/“c-f(x)_7  -  c*Af(x) 


By  rewriting  this  relation  from  right  to  left,  we  obtain: 


cAf(x)  -A/”cf(x)_7* 


This  is  a  mathematical  formulation  of  the  second  part  of  the 
assertion  under  consideration. 

'3.  For  an  entire  function  of  the  kth  degree: 
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y  -  Aq (x  -  xi)k  +  AjCx  -  x^)1*"1  +...+  Akl(x  -  xi)  +  Ak; 


/  ■  th  • 

■/  the  k  differences  are  identical. 

* 

Let  us  confine  ourselves  to  the  case  of  k  *  2,  so  that 


-1-  i 


Aq(x  -  x t-  A^(x 


Xf)  +  Ag ■ 


Taking: 


*i+l  "  xi  +  h; 


Xi+2  "  xi  +  2h; 


xi+3  “  xi  +  3h* 


we  obtain  from  the  formula  for  the  entire  function: 


yi  ”  V 


yi  "  yi4-l  -  Vi  -  A0h  +  Alh: 


yi+i  -  Aoh  +  Aih  +  A2  * 


yi+l  "  yi+2  -  yi+l  “  3A0h  +  AlhJ' 


y i+2  ”  4A0h  ^  2\h  +  A2  * 


Ayi+2  “  yi+3  “  yi+2  “  5Aoh  +  Aih; 


yi+3  "  9A0h  +  -3A1h  +  A2; 


and  finally: 


A2yt  -  Ayi+i  -  AyA  -  2AQh2; 


*  ?i+l  -  *>*1+2  "  Ay 


i+l  "  2A0h 
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The  differences  A~y.  and  A  y.  ,  are  found  to  be  identical 

1  irl 

(because  and  h  are  constant  numbers),  which  is  what  we  set 

out  to  prove. 

■  '  '  4)  Determination  of  Coefficients  of  an  Entire  Function. 

Let  us  undertake  the  task  of  expressing  the  coefficients  of 
ah  entire  function  in  terms  of  differences  of  this  function,  lo 
do  this  for  the  case  of  k  »  2,  we  shall  make  use  of  the  following 
relations: 


+  Ajh; 


2A  h2. 
0 


Substituting  the  value  of  A^  from  the  last  relation, 


namely: 


into  the  equation  forAy^,  we  find: 

a  1  .2 

Ayf  -  ---  A  yA  +  Ajh, 

whence 

A.  -  -1*  (2Ay .  -  A2y  ) 
1  2h  .  1  1 

Thus,  finally: 

.  2^i  - 

Ao  ■  -s-!  A,  ■  - ; - : 

v  2h2  1  2h2 
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The  eoef f icients  . lor  an  entire  third  doBr.ee  function  can  be  found 


in  a  similar  manner: 

A3 

#  A  yi 

Ao  "  — 7  ’ 
u  6h3 


S’ 


Ao  “ 


6Ay.  -  3A“yj  +  2h  y i 


Ao  -  y . . 


These  coefficients  will  be  needed  subsequently  for  deriving  the 
formula  of  the  interpolation  function.  It  is  of  interest  to  note  the 
fact  that  the  coefficients  of  the  entire  function  are  here  expressed 
in  terms  of  differences  with  the  same  subscript  i. 

The  same  method  may  be  used  to  determine  these  coefficients 
in  terms  of  differences  with  different  subscripts.  Thus,  for 
example,  by  taking  tj,  ti_1,  ti2,  ti3  (instead  of  t^  *i+y  ti+2, 
ti+3>*  we  sha^  find  the  following  relations: 


A0  “ 


A  Yi-3 


^  V  +  A3  V  n 

-  >i"2  yi-3. 


.  **1.1  +  3a\-2  +  2A\-3 
2  to  1 


A_  -  y ,  . 


These  relations,  which  the  reader  himself  will  be  able  to  derive 
directly  by  using  the  same  method,  will  find  an  application 
in  the  derivation  of  working  formulas  for  the  numerical  integration 


of  equations. 


5)  The  Practical  Value  of  Differences. 


Differences  are  employed  for  the  following  purposes: 
a)  Determination  of  intermediate  values  of  a  function. 
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b)  Computation  of  intermediate  values  of  the  argument. 

c)  Factual  determination  of  derivatives  of  various  orders  of 
a  function. 

d)  Determination  of  definite  integrals. 

e)  Numerical  integration  of  ordinary  differential. equations. 

The  universal  character  of  thetfe  differences  used"  in  the 

operations  enumerated  above  is  reflected  in  the  fact  that  it  is 
perfectly  immaterial  whether  the  function  has  been  stated  in  the 
fora  of  a  tabulation,  a  diagram,  or  a  formula. 

It  has  already  been  established  that  the  k**V  difference  of  an 
entire  function  of  the  k**1  degree  is  constant.  It  is  not  difficult 
to  show  that,  conversely,  if  the  differences  of  the  k**1  order  of 
a  certain  function  are  constant,  the  latter  is  an  entire  function 
of  the  kth  degree. 

This  proposition  constitutes  the  cornerstone  of  the  utilization 
of  differences  in  all  of  the  operations  indicated  above,  for  the 
following  reason.  If  the  formulation  of  the  table  of  differences  of 
a  certain' function  shows  -that  the  k^k-diffei'ences  are  almost  constant, 
we  have  the  right  to  replace  our  function  by  an  entire  function  of 
the  kth  degree  and  then  subject  the  latter  function  to  all  the 
necessary  operations. 

It  is  this  entire  function  which  is  designated  as  the  interpolation 
function.  A  general  expression  for  it  must  be  derived.  When  k  -  3, 
we  have : 

o  ~  ~  2  ~ 

y  •  Aq(x  -  x^)3  +  -Aj/<x  -  xi).  +  AgU  -  xA>  +  A3. 
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The  interpolation  function  is' -thus-found  to  be  expressed  in 
terms  of  differences  of  the  given  function  up  to  and  including  those 
of  the  third  order,  with  all  these  differences  having  one  and  the 
same  subscript. 

% 

In  replacing  any  function  (regardless  of  the  manner  in  which  it 
is  stated)  by  this  interpolation  function,  it  -is  categorically  . 
F-TS-7327-RE  625 

" .  -r  r-‘  :  »  ■■ 

rs  *  t  *  t 


CHARLES  A.  MEYER  &  CO.  Ini 


imperative  to  dir.ect  attention,  on  the  basis  of  the  table  of 

*  * 

differences  of  the  given  function,  to  the  character  of  the  variation 
of  these  differences;  it  is  necessary,'  and  triply  necessary,  that 

the  differences  be  nearly  constant. 

.  •  « 

Obviously  this  condition  will  bo  satisfied  the  better  the  smaller 

the  absolute  values  of  the  differences  of  any  order. 

Rule.  In  order  to  reduce  the  differences  of  any  order,  the • 

step  must  be-  reduced. 

/ 

f  This  rule  results  from  the  following  relation: 

Aky  -  y(k)  •  hk  +  h*  •  £-  hk(y^  +  t). 


where  y^  is  the  k^-order  derivative  of  x,  and  t  is  an  infinitely 
small  quantity  of  the  highest  order.  The  relation  itself,  which 
is  derived  at  the  proper  stage  of  the  differential  computation, 
represents  a  generalization  of  the  better-known  relation  between 
the  increment  of  a  function  and  its  differential: 


Ay  «  dy  +  ht 


Ay  -  y’h  +  h€  «  h(y'  +  E). 


Thus,  as  the  step  is  halved,  the  first  differences  decrease 
approximately  by  a  factor  of  two,  the  second  differences  decrease 

by  a  factor  of  four,  and  the  third  differences  decrease  by  a  factor 

.  / 

V  *  -  -  -  

of  eight. 
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allies  of  a  Given  Function. 


To  determine  the  intermediate  values  of  a  given  function, 
regardless  of  how  it  is  stated,  it  is  necessary  to  formulate  a 
table  of  differences  of  this  function  and,  having  made  certain 
that  the  second  differences  are  constant  or  are  very  small  in 
comparison  with  the  particular  values  of  the  function  itself,  replace 
it  by  an  entire  function  of  the  second  degree,  i.e.,  by  an  interpolation 
function  (105) : 


y-yi+^  •Ayi-~§(1-  5)A“yi. 


(106) 


where  y^  designates  the  -tabular  value  of  the  given  function 
corresponding  to  the  next  smaller  value  of  the  argument  nearest 
to  that  value  of  x  for  which  the  intermediate  value  of  the  function 
is  sought. 

Knowing  the  intermediate  value  of  x  for  which  the  intermediate 

value  of  the  given  function  is-  sought,  and  having  obtained  from  the 

table  of  differences  of  this  function  the  nearest  value  of  x  ,  we 

\  * 

can  determine  the  value  of  the  coefficient  of  interpolation  £  and 
the  unknown  value  of  y  by  means  of  formula  (106). 

Example  1.  Given  the  following  table  of  pressure  impulses: 


where  the  impulses  are  expressed  in  kg^dm"  -sec,  determine  the 
pressure  impulse  for  t  -  0.020556. 

Solution.  Lot  us  formulate  the  table  of  differences: 

Table  4-d. 


t-103 

20 

21 

22 

1 

71.0 

80.6 

to 

to 

• 

A I 

9.6 

11.5 

a2i 

1.9‘ 

Let  us  find  the  coefficient  of  interpolation: 


t  -  t4  0.020556  -  0.020 

5-  - -  -  -  -  0.556. 

h  0.001 


The  desired  intermediate  value  of  the  function  is: 

I  -  I.  -  4*(1  -  £)A2I.  -  71.0  +  0.556  *  9.6  —0.556(1 

2  1  •  2 


-  0.556)1.9, 


since 

$-  0.556;  IA  -  71.0;  -  9.6  and  -  1.9. 

We  finally  obtain: 

I  «■  71.0  -f  5.33  -  0.25  »  76  kg,dm_2‘sec 
,  7)  Computation  of  Intermediate  Values  of  the  Argument. 

Cases  are  sometimes  encountered  in  practice  of  a  contrary  nature 
in  which  an  intermediate  value  y  of  the  function  is  given  and  it  is 
required  to  find  the  coefficient  of  interpolation  £  and  the  intermediate 
value  of  the  argument  x.  Such  an  operation  is  sometimes  called  an 
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inverse  interpolation. 

Then,  discarding  the  last  term  in  formula  (106),  we  obtain  as 
a  first  approximation: 


y  -  .y. 

,  s. - 

y  1  AyA 

As  a  second  approximation,  we  find  for  ?  from  the  same  formula 


(106) 


y  -  yi 


Ayi  "  T(1  " 


Let  us  assume  in  the  right-hand  side  of  the  above  expression 


%  -  5^,  then  in  the  second  approximation: 


:*r- 


y  -  yn 


Ay,  -  -|(1  -  S1)A2y1 


where  5^  is  already  known  from  the  first  approximation. 

Sometimes  the  following  expression  is  used  as  the  third 


approximation : 


*3- 


y  -  y. 


Ay,  ,'L(1  -  UA,  +  |(1  -  y<2  -  ?,)A3y 
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This  formula  can  bo  easily  obtained  from  the  general  expression 
for  the  interpolation  function  (105)  by  placing  £*outside  the 
parentheses,  determining  a'nd  substituting  in  the-  right-hand  side 

|  ■  %  known  from  the  second  approximation. 

, _  ~  . 

.  The  desired  value  of  x  will  be: 


x  -  x .  +  \  h  or  x  -  x  +  % oh. 

1  m  1  O 


Example  2.  We  have  the  following  table  for  the  relative 
portion  vy  of  the  charge: 

Table  4-e. 


t-103 

20 

21 

22 

Y  . 

0.034 

0.042 

0.054 

It  is  desired  to  find  the  value'  of  t  corresponding  to  the  inflow 
of  gases  ^  -  0.038. 

Solution.  Let  us  formulate  the  following  table  of  differences: 

Table  4-f. 


t-103 

20 

21 

22 

•  ^ 

A 

A2h» 

0.034 

0.008 

0.004 

0.042 

0.012 

0.054 

We  obtain  as  a  first  approximation: 

*  • 

^  -  w.  0.038  -  0.034  i 

w  -  _ — i  -  — - -. 

•  AH'i  0.008  2 

We  find  as  a  second  approximation: 
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M*-  H». 


JL 


AT;  *  -^d  -  51)A%1 


0.038  -  0.034 
0.008  -  i.  -i-0.004 


--  -  0.556. 
7 


The  desired  value  of  the  argument  is: 


t  -  t  +  5  h  -  0.020  +  0.556(0.024  -  0.023)  -  0.020556. 

.  *  u 

8)  Numerical  Differentiation  of  Functions. 

For  this  operation  it  is  necessary,  first  of  all,  to  formulate 
a  table  of  differences  of  the  given  function  and  to  make  absolutely 
certain  that  the  differences  of  the  third  order  are  nearly  constant. 
This  fact  makes  it  possible  to  replace  the  given  function  (regardless 
of  the  manner  in  which  it  is  stated)  by  the  interpolation  function 
(105) :  . 

y  -  y  +  5  •  Ay.  -  —5(1  -  5  )A2y  +  —  5(1  -  5)  (2  -  S)A3y.. 

1  2  1  6  1 

From  this  we  find  the  desired  derivative: 


dy 

dx 


d?  1  d5 

-  Ay ,  +  —  25 - 

dx  i  2  \  dx 


dj_ 

dx 


Asyi  +  —  1 2  ~ 
1  6  \  dx 


65  —  + 
dx 


+  35- 


d5_ 

dx 


A3y 


i' 


n  O 

keeping  in  mind  that  the  differences  y^,  Ay^,  A  y^,  and  A°y^  are 


certain  constant  numbers. 
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If  the  step  h  is  not  a  prime  number  (for  example, 


constant  h  may  be  introduced  under,  the  sign  of  the  di 
any  desired  order.  Therefore,  upon  introducing  into 
the  auxiliary  function: 


— ),  the 
738 

fference  of 


the  process 


*  -  -fy* 

h 


we  obtain 


A.  1  A 2  12  .3  13 

Aq>  -  — Ay .  :  A  4> .  -  . —  Ay.;  A  $  *  —A  y . , 
i  h  *  i  h  1  i  h  i 

and  the  working  formula  for  numerical  differentiation  acquires  its 
final  form: 

( 

l  dx 


1  ,2.  1  3, 

A4.  -  —  A  i.  +  — A  $>  , 

1  2  1  3  i 


Example  3.  We  have  the  following  table  for<y,  the  relative  portion 
of  the  charge: 


Tab 


e  4-g. 


t.io3 

21 

21.5 

22.0 

22.5 

*103 

42 

47 

54 

64 

diy 


It  is  desired  to  find  the  rate  of  gas  formation  —  at  the  instant 

dt 


t  -  0.0210. 

Solution.  Let  us  formulate  the  following  table  of  differences: 
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Tabic  4-h, 


/ 

y 


t*  103 

21 

21.5 

22.0 

22.5 

H'yio3  • 

42 

47 

54 

64 

Avy  •  io3 

5 

7 

10 

A  2^. 103 

2 

3 

Zt^.jQ3 

1 

3 

Since  the  given  values  are  t^*10  «  21  andvy^‘10  ■■  42,  it 

follows  that: 

A  (>|»  •  103)  -  5;  A2(H’i,103)  -  2;  A3('Vi*  103)  -  1. 


\ 

The  first  of  the  working  formulas  for  the  derivative 


4L. 

dx 


will  make  it  possible  to  determine  the  rate- of  gas  formation  (h  - 


0. 5* 10°) : 


d(fl03)  .  1 

dt  h 


A^-IO3)  -  ~A2('Vi*103)  +  ~A3(4'j  *  103) 


from  which: 


5-1 

0.5 


"  8 


sec 


t-0.021 

9)  Computation  of  Definite  Integrals. 


F-TS-7327-RE 


634 


CHARLES  A!  MEYER  &  CO.  Inc. 

_ _  an  ate  Aaa 


y 


Fig.  152 


If  it  is 


/ 


-  Determination  of 


I 


ydx  as  a 


function  of  x  from  the 


curve  y  -  f (x) . 

necessary  to  find  any  definite  integral: 


a 


<* 


where  the  limits  of  integration  a  and  b  are  finite,  this  operation 
is  equivalent  to  computing  the  shaded  area  shown  in  fig.  152,  which 
is  limited  at  the  top  by  a  curve  representing  the  function  whose 
integral  is  to  be  found.  To  compute  this  area  Y,  the  usual  method 
is  applied  first:  the  interval  of  integration  (b  -  a)  must  be 
divided  into  n  equal  portions,  lines  normal  to  the  0-x  axis  are 
then  erected  at  the  po.ints  of  division,  and  the  unknown  area  is 
divided  into  elementary  areas: 


AY 


O’ 


AY, 


AY 


i' 


AY 


n-1 


The  problem  is  then  reduced  to  the  determination  of  these 
elementary  areas.  For  this  purpose  it  is  sufficient  to  replace  the 
portion  of  the  curve^-f  (x^)  corresponding  to  the  elementary  area 
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from  which,  after  substituting  the  extreme  values  of  we  obtain: 

*'i ->(>.■•  4  *’i-b**'i*  a  *■>.)• 


This  is  the  working  formula  employed  in  finding  a  definite 

/ 

integral  by  the  method  of  numerical  integration. 

To  determine  the  unknown  integral,  it  remains  necessary 
to  summate  gradually  and  successively  the  individual  elementary 
areas : 


Yx  -  Y0  +AYq; 


Yi  "  Yi-1  +AYi-r 


Y2  -  Yx  +AYX; 


Y  -  Y„  ,  + 
n  n-i  n-1 


Example  4.  We  have  the  following  table  for  a  function  to  be 
integrated;  #  | 

_ Table  4-i. _ 

flO3  8  12  16  20 

p,  kg'cra“2  26  32  48  88 

It  is  desired  to  find  the  elementary  area  Alj,  the  subscript  i  being 
applied  to  the  first  column.  \ 


Solution.  Let  us  formulate  the  following  table  of  differences: 

Table  W j  . _ 


t  •  103 

\ 

8 

12 

16 

p 

Ap 

~2S~ 

«\ 

16 

~W5 

40 

F-TS-7327-RE 


CHARLES  A.  WIVES  &  CO.  tru  . 


Table  4-j  (Cont'd.) 


a2p 

10 

24 

- 

- 

aV 

14  . 

- 

- 

Since : 


h  -  4- 10“3 ; 


yi  -  26;  A yi  -  6;  A2yA  -  10;  A3yi  -  14, 


the  elementary  area  under  consideration  will  be: 


AYi  - h  [yi +  iAyi  -  +  ^3yi 

-  4-10“3  [26  +  3-  —  4-“]^  4-10-3(26 
^  12  24  / 


+  3  -  1  +  1)~ 


~  4-10~3*29  ~  0. 116  kg-sec*cm“2  ss  11.6  kg- sec -dm""2. 


10)  Determination  of  .Pressure  Impulse  from  a  Pressure- 

Bomb  Test. 

As  is  already  known  from  the  preceding  course  in  internal 

ballistics,  the  pressure  impulse  is  expressed  by  the  following  relation: 

t 


-  j  pdt; 


It  is  clear  that  in  order  to  determine  it  from  fi  'bomfPtest 
yielding  directly  the  (p,  t)  curve,  it  is  necessary  to  compute  a 
definite  integral,  for  which  the  pressure  of  the  powder  gases  p  is 
the  function  being  integrated  and  the  time  t  is  the  independent 
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variable.  This  can  be  accomplished  of  course  by  means  of  one  of 

4  *  *# 

the  usual  methods  of  quadratures.  However,  it  is  usually 
necessary  to  have  a  durve  for  the  pressure  impulse  as  a  function 
of  the  gas  inflow^;  in  other  words,  in  addition  to  the  pressure 
impulse  I ^  at  the  end  of  the  burning  of  the  powder,  it  also 
y  becomes  necessary  to  find  a  series  cf  intermediate  values  for 

this  quantity,  which  will  correspond  to  intermediate  values  of  the 
gas  inflow.  This  problem  is  solved  most  simply  by  the  method  of 
numerical  integration. 

Example  5.  Given  the  pressure  curve  (p,  t)  represented  by 
Table  4,  find  the  correlation  (I,  y)  by  means  of  the  table  correlating  y 


with  t. 


_ Table  4-k  -  .Experimental  Table  of  as  a  Function  of  t. _ _ 

t-103  0  4  8  12  16  18  20  21  22  23  23.5  24  24.5 


k>-10J  0  2 


6  14  22  34  42  54  79  101  127  155 


t-103  25  25.5  26.0  26.5  27.0  27.5  28.0  28.5  29.0  29.5  30.0  30.5  31.0 

4>-103  1  86  221-  2  6  0  304  3  54  409  4  70  53  9  619  711  813  9  1  9  1  000 


Solution.  We  perform  all  computations  on  the  working  form  in 
Table  5,  without  considering  third  differences,  and  designating: 

•  „  1  12  1  1  2 

y  +  — Ay  -  — A  y  -  p  +  — Ap  -  —  A  p: 

2  12  2  12  y 


from  which  it  follows  that: 


M.  -  AY.  -  h  y.  +  —Aya  -  ~A2y.  ]  -  h-£. 
x  i  ,  *  9.  1  19  i< 
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Table  5  -  Computation  of  Pressure  Impulse 


j  j*Mweicw^»fga.yn3ttSBfi*£y  - — 


x 


«r 


*w 


The  computations  in  Table  5  are  conducted  in  rows,  proceeding 
from  left  to  right  in  each  row  and  downward  from  one  row  to  the 
next,  in  the  following  manner. 

3 

1)  The  first  row  is  filled  with  values  of  the  argument,  t*10  . 

2)  The  numbers  for  the  pressure  are  taken  from  Table  4  (page 

616).  - 

3)  From  each  number  in  the  second  row  there  is  subtracted 
the  preceding  number  in  the  same  row,  and  the  result  is  written 
under  the  left-hand  number  of  this  pair: 

23  -  21  -  2;  26  -  23  -  3;  etc. 

4)  The  numbers  in  the' fourth  row  are  obtained  in  the  same 
manner  as  in  the  preceding  row: 

3-2«  1;  6-3-3;  etc. 

5)  The  numbers  of  the  second  row  are  repeated. 

6)  The  numbers  in  the  third  row  are  halved. 

7)  The  numbers  in  the  fourth  row  are  divided  by  12  and  written 
with  the  opposite  sign. 

8) .  The^ numbers  in  the  three,  preceding  rows  are  added. 

9)  The  \numbers  in  the  . eighth  row  are  multiplied  by  the  step  h, 


the  decimal  point  being  correctly  placed  (cf.  Example  4)  to  convert 


o  \  _o 

kg  .sec*cm“'<!  into  kg'sec’dm**  . 

10)  To  obtain  the  next  succeeding  value  of  I^,-  it  is  necessary 
to  add  to  the  preceding  value  of  the  corresponding  increment 

Mi_l: 


Ii  ”  xi-l  +  AIi-l» 
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or,  in  other  words,  it  is  necessary  tc  add  the  two  numbers  in  the 
ninth  and  tenth  rows  of  the  preceding  column: 

0  8.8  -  8.8; 

8.8  -f  10.0  -  18.8,  etc. 

By  taking  the  values  for  the  gas  inflow  ^  and  the  pressure 
impulse  1  corresponding  to  the  same  instants  t,  it  is  easy  to 
obtain  the  desired  correlation  between  1  and  (cf .  Table  6). 


Table  6  -  Correlation  between  I  and  H*  for  Test. 


o 

I  kg*sec*cJm 

0 

9 

19 

30 

45 

56 

71 

81 

92 

107 

117 

129 

144 

162 

v^-103  . 

0 

2 

3 

6 

14 

22 

34 

42 

54 

79 

101 

127 

155 

186 

-2 

I  kg*sec*dm 

184 

209 

313 

• 

359 

412 

473 

544 

625 

718 

770 

823 

3 

*  *10 

221 

260 

409 

470 

539 

613 

711 

813 

919 

967 

1000 

Table  6  will  be  needed  subsequently  for  the  solution  of  the 
principal  problem  of  internal  ballistics  (determination  of  the  curves 
for  the  speed  of  the  projectile  and  for  the  pressure  of  the  powder 
gases  as  a  function  of  the  path  of  the  projectile).  This  solution 
must  be  prepared  by  a  discussion  of  the  necessary  theory,'  which  will 
be  undertaken  in  the  next  section. 

•  2.  NUMERICAL  INTEGRATION  OF  DIFFERENTIAL  EQUATIONS 
11)  Numerical  Integration  of  First-Order  Equation 
We  have  such  an  equation  in  its  general  form: 

F(y*  y',  x)  -  o. 
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It  is  composed  of  the  following  elements:  the  independent 
•  >, 

vari:«ble  x,  the  unknown  function  y,  and  its  fjrst  derivative  y'. 
•It  is  necessary,  first  of  all,  to  determine  this  derivative  from 
the  general  equation: 


y'  -  f(y,  x); 


It  is  this  equation  which  is  solved  by  the  method  of  numerical 
integration,  regardless  of  the  form  of  the  function  which  constitutes 
its  right-hand  side.  In  practice,  what  is  usually  found  is  not 
the  total  integral  of  the  last  equation  (in  which  case,  of  course, 
numerical  integration  cannot  be  employed),  but,  rather,  a  partial 
integral,  in  consequence  of  which  it  is  necessary  to  indicate  the 
’’initial  conditions",  i.e.,  the  values  of  y,v  and  x  .  In  addition,' 
the  value  of  x^  at  the  end  of  the  finite  region  of  integration  must 
also  be  known.  Under  these  conditions  the  desired  function  will  be: 


or 


y 


f(y,  x)dx. 


f 

! 

t 

1 


The  problem  is  consequently  reduced  to  the  finding  of  a  definite 
integral  whose  characteristic  feature  is  such  that  the  part  of  the 
function  to  be  integrated  is  played  by  the  derivative  of  the  function. 
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•It  would  seem  that,  as  on  the*  previous  occasion  involving  the  finding 
of  a  definite  integral,  it  might  be  possible  to  perform  computations 

by  rows  and  to  make  use  of  a  working  formula  containing  differences 

/ 

of  an  auxiliary  function  with  they4:.me  subscript  i.  In  actual 
practice,,  however,  this  procedure  cannot  be  adopted.  As  a  matter 
of  fact,  the  last  formula  shows  that  the  function  to  be  integrated 
includes  the  original  function  itself,  and  we  ttyus  obtain  a  "vicious 
circle;" -in  order  to  determine  the  particular  value  y^  of  the  function 
we  must  know  the  particular  value  4^,  and  in  c/rder  to  determine 
this  particular  value  of  the  auxiliary  function: 

/ 

*i  -  f(yi>  xi)h 

we  must  have  the  value  of  y^  which  is  a  constituent  of  the  auxiliary 
function. 

It  is  now  clear  that  in  performing  numerical  integration  of  the 
equations  it  is  not  possible  to  perform  computations  by  rows,  it 
being  necessary  instead  to  proceed  gradually,  step  by  step,  performing 
the  required  computations  downward  along  each  column  and  to  the 
right  from  one  column  to  the  next. 

Let  us  further  assume  that  the  values  4*.  .  0,  4>,  and  4>. 

of  the  auxiliary  function  have  already  been  found.  With  their  aid 
it  is  possible  to  compute  the  following  differences. 
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Table  6-a. 


•  $ 

1-2 


A4>  U* 


H^i-2 


A2*  A  -^..3  A2^_2 
A'\  ^3^i-3 


The  above  table  shows  that  in  order  to  determine  Ay^,  it  is  not 
possible  to  make  use  cf  the  following  formula  for  determining  definite 


integrals: 


-*i +  -Wi  -  +  nAV 


2  ^ 

because  the  differences  A  4*^  and  A  4*^  are  not  yet  contained 

in  the  table.  Consequently,  it  becomes  necessary  to  derive  an 
additional  formula  containing  the  differences: 

A*i-1’  :A\-2’  aVs' 


already  contained  in' the  table. 

For  this  purpose,  as  on  the  previous  occasion,  we  shall  replace 
the  function  to  be  integrated,  y'  -  f(y,  x),  by  the  following 
interpolation  function: 

3  2 

y «  A  (x  -  x^)  +  Ax(x  -  xA>  +  A2(x  -  xA)  +  A^v 


and  then: 
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f(y,  x) dx 


J  rv* 


c  -  Xi)‘  +  A  j  ( x  -  xi)“  +  A2(x  -  xj  +  Ao  7dx  - 


7ao(x 


xi)‘1  +  -—AjCx  -  xi)‘J.+  •—  Ag Cx  -  xA)2  +  A3(x  -  xA) 


whence 


1  .  u4  1  .  u3  .  1  »  u2  ,  .  . 

4  A0h  +  3  Alh  +  2A2h  A3h> 


because 


*i+i  -  *i  -  »• 


In, contrast  with  the  preceding  case,  we  shall  substitute  here 


the  following  values  for  the  coefficients: 


a 

0  6h3 


A1  - 


A2y!_2  +  A3y^3 


A2  “ 


6Ay.*_1  +  3A2y^_1  +  2A3y!_3 


f  -  y.’* 


which  are  expressed  in  terms  of  differences  of  the  function  y  uuder 
the  integral  sign  provided  with  the  required  subscripts. 

Following  obvious  transformations,  we  obtain: 
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It  is  usually  convenient  to  make  use  of  the  following  auxiliary 
function: 

$  -  hy1, 

the  derivative  y'  being  computed  in  accordance  with  the  given  equation: 

y*  -  f  (x,  y). 


y: 


Then: 

hy’-*!;  hAy’,j  **\-2 

and  the  working  formula. for  the  numerical  integration  of  ordinary 
differential  equations  of  the  first  order  will  assume  the  following 
form: 

AX!  -  *i  +  +  ^A2ii.2  +  -§-A3*i-3; 

This  formula  now  includes  differences  of  the  auxiliary  function  $ 
provided  with  precisely  those  subscripts  which  are  required  in 
•  accordance  with  the  table  of  difference  presented  above. 

It  will  be  seen  from  the  same  table  that  in  order  to  determine 
v  Ay^,  it  becomes  necessary  to  take  differences  from  different  columns, 
which,  of  course,  is  inconvenient.  For  this "reason  it  is  desirable 
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to  modify  the  system  of  writing  the  differences  in  such  a  manner  as 
to  have  in  the  same  column  differences  with  different  subscripts, 
na  me  1  y ; 


Table  6-b. 


4>  ■ 

+  .  o 

1-3 

•  4> 

j-2 

<K 

l 

A* 

■ 

A$.  „ 
1-3 

Aii-2 

A4>.  . 
l-l 

A2* 

*2*i-3 

A2*  1-2 

A34> 

In  such  a  procedure  each  column  will  contain  those  differences 
which  are  necessary  for  the  application  of  the  new  working  formula 
for  numerical  integration. 

Rule  for  writing  new  differences.  From  the  right-hand  number 
in  a  row  there  is  subtracted  the  left-hand  number  in  the  same  row, 
but  the  result  is  written  under  the  right-hand  number. 

Finally,  from  the  last  table  of  differences  there  follows  the 
.most  disagreeable  characteristic  feature  of  the  numerical  integration 
of  equations:  in  order  to  determine  Syx,  it  is  necessary  to  know  the 
differences : 

A4»  .  i ,  A  4s .  q  ,  a  4>  •  nt 

1— X  1“*  1— J 

i.e.,  to  have  the  following  particular  values  of  the  function: 

4>  ,  4> ,  9,  4> 0. 

^  i-1  X~£>  i-o 

These  will  become  known  in  the  course  of  the  process;  but  at 
the  start  of  the  integration  only  one  particular  value  of  this  function, 
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namely^,  is  Known  from  the  originally  stated  equation  for  y* 
In  order  to  find-  -• 


$  ~l>  $_2*  *  o 


it  is  necessary  to  employ  special  supplementary  methods,  of  which 
the  .most  commonly  used  is  the  method  cf  successive  approximations. 
The  essence  of  this  method  resides  in  the  fact  that,  ?it  the  start 


of  the  computations,  there  is  adopted  a  gradual  and  stepwise  advance, 

"  >■» , 

each  new  approximation  permitting  one  additional  step,  in  which 

those  differences  that  have  already  appeared  earlier  are  utilized. 

y 

There  exist  several  variants  cf  this  method,  one  of  which,  the 
most  exact  one,  can  be  best  studied  by  the  aid  of  a  concrete  example. 

Example  6.  Solve  by  the  method  of  numerical  integration  the 
equation: 


y  -  x  -  y 

over  the  interval  0;  0*5,  with~a  step  of  h  «  0.1,  if  y  -  1  when 


x  -  0. 


Solution.  All  the  computations  are  presented  in  Table  7. 
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iX&ratfrn  approximation  Third  approximation  Normal  computations 


First  Approximation 


First  Co Uimn 

13)  This  is  the  leading  row  and  is  filled  first;  in  the  given 
case  yy  m  1.0000  from  the  conditions  of  the  problem. 

‘  1)  Xq  -  0  from  the  conditions  of  the  problem. 

2)  -yg  -  -1.0000  in  accordance  with  the  thirteenth  row. 

3)  yg  -  -1;  in  accordance  with  the  given  equation  for  y', 
we  add  the  numbers  in  the  first  and  second  rows  above. 

4)  -  -0.1000;  the  number  in  the  third  row  is  multiplied' 
by  the  step  h  «  0.1,  since: 


*0  *  hJo 

The  fifth,  sixth  and  seventh  rows  are  not  filled  in  for  lack 
of  data.  Then: 

8)  «*  -0.1000;  the  number  in  the  fourth  row  is  repeated. 

The  ninth,  tenth  and  eleventh  rows  are  not  filled  in  for  lack 
of  data.  Thereupon: 

12)  A  }’q  -  -0.1000;  in  this  case  it  is  necessary  to  repeat 
the  number  in  the  eighth  row. 

Second  Co lumn  (of  the  First  Approximation ) 

13)  Leading  row;  y^  -  0.9000;  we  add  (algebraically)  the  numbers 
in  the  twelfth  and  thirteenth  rows  of  the  first  column,  since: 

y-i  -  y0  +  *y0; 


1)  -  0.1;  we  add  the  step  h. 


F-TS-7327-RE 


651 


CHAr.Li  S  A.  .V.LYER  &  CO.  Inc. 


El 


ft- 


iW '"— ' 

'JI 


1 


.2)  -y^  -  -0.9000  in  accordance  with  the  thirteenth  row. 

#  % 

3)  y^  «  -0.8000;  in  accordance  with  the  given  equation  for  y' 

we  add  the  numbers  in  the  first  and  second  rows  of  this  column  (above). 

4)  *  -  -0.0800;  the  number  in  the  third  row  is  multiplied 

by  the  step  h  -  0.1. 

Thus  a  number  appears  in  the  fifth  row;  A4*q  -  200,  for  which 
we  subtract  from  the  number  (-0.0800)  in  the  fourth  row  above  the 
number  (-1.0000)  at  the  left  in  the  same  row,  omitting  the  zeros 
at  the  left  for  greater  convenience. 

Vedge  of  First  Approximation. 

It  consists  merely  of  a  single  number  and  is  surrounded  by 

a  heavy  line.  .  To  obtain  it,  we  postulate  that 


“  200. 


Second  Approximation 


First  Column. 


13)  Leading  row:  $  -  1;  1)  x  -  0.  ' 

V  U  • 

The  second  and  third  rows  are  not  filled,  as  this  is  no  longer 
necessary  (cf.  first  column  of  first  approximation).  Ihen: 

4)  -  -0.1000,  as  in  the  same  row  of  the  first  approximation. 

-  Omitting  the  fifth,  sixth,  and  seventh  rows,  we  write: 


8)  -  -0.1000;  the  number  in  the  fourth  row  is  repeated. 


9) 


1 


■A$  “  100;  we  take  the  difference  A<f>  from  the  wedge 

2  “1 •  -1 

of  the  first  approximation. 

The  tenth  and  eleventh  rows  are  not  filled  for  lack  of  necessary 

data. 
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12)  Ayy  -  -O.ODOU;  we  add  the  numbers  in  the  eighth  and  ninth 
rows  above. 

The  second  column  of  this  approximation  is  formed  in  an 

/ 

analogous  ma i/ner,  with  the  only  difference  that  numbers  appear 
in  the  first  and  second  rows  (cf.  the  number  in  the  thirteenth  row 
of  this  column)  and  in  the  fifth  row,  for  which  we  subtract  from 
the  number  (-0.0810)  in  t'li£  fourth  row  the  number  (-0.1000)  in  the 
same  row,  but  on  the  left,-  The  number  m  the  leading  row  is  obtained 
by  adding  the  numbers  in/the  twelfth  and  thirteenth  rows  of  the 

■"s.. 

preceding  column. 

The  third  column  is  distinguished  from  second  in  that  a  number 

/ 

(-17)  appears  in  the  sixth  row,  for  which  we  subtract  from  the  number 
173  in  the  fifth  row  the  number  190  at  the  left,  and  another  number 
(-8)  appears  in  the  tenth  row,  for  which  the  number  in  the  sixth 
row  is  multiplied  by 


5  io 
12  24* 

Filling  in  the  Wedge  of  the  Second  Appro xima t i o n . 


In  the  sixth  line  we  repeat  twice  the  number  (-18),  noting  that: 

A2$  «  A2$  •  -  a24>  «  -18. 

-2  -1  0 


The  number  208  in  the  fifth  row  of  the  first  column  of  the  second 
approximation  will  appear  in  accordance  with  the  definition  of  the 
second  differences: 

-A* 0  -**-1’ 

from  which:  / 
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A*  -  A4>_0  -  A24,1  -  li»0  -  (-18)  -  208. 

Third  Approximation. 

In  the  first  column,  after  the  thirteenth  and  first  rows  are 

filled,  the  places  in  the  second,  third,  fifth,  sixth,  and  seventh 

rows  are  left  blank.  Then; 

8)  The  number  in  the  fourth  row  is  repeated. 

0 )  —  A$  ■*  104;  the  number  208  in  the  wedge  of  the  second 

2  —I 

approximation  is  halved. 

5  2 

10)  — A  ♦  "  -8;  the  number  (-18)  in  the  wedge  of  the  second 

12 

approximation  Is  multiplied  by 

24 

The  eleventh  row  is  omitted.  Thereupon: 

12)  The  numbers  in  the  eighth,  ninth,  and  tenth  rows  above  are 

added : 


-0.1000  +  0.0104  -  0.0008  -  -0.0‘j04. 

In  the  second  column,  after  the  thirteenth,  first,  second, 
third,  fourth  and  fifth  rows  are  filled,  the  sixth  and  seventh  rows 
are  omitted.  Then: 

8)  The  number  in  the  fourth  row  is  repeated. 

9)  —  **  y5;  the  number  in  the  fifth  row  of  tbis  column 

2  0 

is  halwed. 

W)  2  •  -8;  the  number  (-18)  in  the  second  column  of  the 

preceding  wedge  is  eultiplied  by 

The  eleweath  row  is  likewise  omitted.  Thereupon: 
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12)  The  numbers  in  the  eighth,  ninth,  and  tenth  rows  above  are 


added . 

In  the  third  column,  only  the  thirteenth,  first,  second,  third, 
fourth,  fifth  and  sixth  rows  are  filled.  There  appears  for  the 
first  time  a  number  in  the  seventh  row: 

-18  -  (-17)  -  -1. 

/ 

Filling  in  the  Wedge  of  the  Ih ird  Approximation. 

We  make: 

*3*-3  "  a3*-2  "  A3*o  ”  -1* 

The  number  (-16)  in  the  sixth  row  of  the  second  column  of  the 
wedge  will  appear  from  the  definition  of  the  third  differences: 

a2*-1  -  -17  -(-1)  ~  -16. 

In  an  analogous  manner,  the  blanks  in  the  rows  of  the  first 
column  of  the  wedge  will  be  filled  as  follows: 

“A2*.!  -aV2  *  ~16  ‘  (_1>  "  "I5; 

2  2, 

-  A*0  -  A  *  x  -  lyO-  (-16)  -  206. 

For  convenience  of  computation,  the  numbers  of  this  wedge  are 
transferred  to  their  places  in  the  columns  of  the  subsequent  normal 
computations.  The  approximations  are  considered  to  be  completed 
because  of  the  cioseness  of  the  results  of  the  third  and  second 
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approximations,  respectively. 


Normal  Computations. 

First  Co ltmn . 

13)  As  before,  y^  »  1.0000;  2)  -  0, 

The  second  and  third  rows  need  not  be  filled.  The  numbers  in  the 
fifth,  sixth,  and  seventh  rows  are  already  in  place.  Thereupon: 

8)  $0  -  -0.  1000;  the  number  in  the  fourth  row  is  repeated. 

9)  — A%  ,  •  103;  the  number  203  in  the  fifth  row  is  halved. 

2  -1 

10)  — a2*  „  «  -6;  the  number  (-15)  in  the  sixth  row  is 

12  in 

multiplied  by  — . 

24 

11)  — A  *  «  •  0;  the  number  (-1)  in  the  seventh  row  is 

8  "*> 

multiplied  by 

12)  Ay^  -  -0.0903;  the  numbers  in  the  eighth,  ninth,  tenth,  and 
eleventh  rows  are  added. 

The  subsequent  columns  are  filled  in  exactly  .he  same  manner. 
Attention  is  once  again  directed  to  the  crder  of  computation  in 
each  column  of  normal  computations. 

a)  First  of  all,  the  leading  thirteenth  row  is  filled  by  adding 
the  numbers  in  the  twelfth  and  thirteenth  rows  of  the  preceding 
column,  because 

-  *i-l  +  A^i-1* 

b)  Thereupon,  the  spaces  in  each  column  are  filled  from  top 
to  bottom  without  omissions. 
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12)  Use  ot  Numerical  I 'ilepnt  t  ion  ol  the  First-Order  Kqu.i  t  ion 

*i_th  Argument  v. 

The  method  of  numerical  integration  must  be  aj.pl  led  to  the 
solution  of.,  the  follo»ing  equation: 


-<L  +  />» 


np  _*! 

v2 

vnp 


where : 


4  -  !ft  -  laif. 


In  the  physical  law  of  pov.de  r  burning  the  correlation  between  vp 
and  v  is  given  by  a  table,  whereas  in  the  geometric  law: 


v+  V«. 

K1k  *  *k 


Since  it  is  also  necessary  as  a  rule  to  find  the  pressure  curve, 
it  is  preferable  to  make  use  of  the  following  equation: 

d  v 

<P»v  ~  -  PS, 


whence 
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Cf  l  i *  fv  •»  .  l  .  ^  U  w  U.  >1  W%- 


where 


P 


8  +  (  ‘ 


It  is  not  difficult  to  see  that  the  auxiliary  function  4>  is 
found  in  the  given  case  from  the  following  relations: 


*- 


fw 


.2 


_ 2£.;  *  .  JL. 

8  1^+1  8  P 


<*) 


As  for  the  working  foraula  of  numerical  integration,  it  has  the  usual 
fora: 


■.4  1  5  2  33 

A,i  -*i  -  +  nA  *i-2  * 


The  purpose  of  the  preliainary  computations  is  to  determine  all  the 
constants: 


i>  i>u 


ts  - 


32.  3)  .*  .  2  I  “ 

g.’  3)  V1p  2  8  « 


t  m 


s 


T  )!  *> 


_1 

6 


_1 

A 


4 


«)  - 


7)  — 


f  1 

=-  +  a  -  T- 


fu 

s 
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lo  determine  the  step  of  integration  h,  the  speed  of  the  projectile 
v j,-  at  the  end  of  powder  burning  must  be  known,  and  then: 


where  n  -  the  number  of  sections  -  is  taken  in  the  range  of  10-40, 
depending  upon  the  density  of  loading  (as  this  density  increases, 
the  number  n  must  also  be  increased). 

While  determining  {  and  p,  the  same  working  fora  may  be  used 
to  find  the  time  of  motion  of  the  projectile  in  accordance  with  the 
following  relation: 

dt  ^  y m  1 
dv  "  8  P  * 


in  the  form  of  a  definite  integral. 

Example  7.  Determine  the  projectile  velocity  and  gas  pressure 
curves  for  a  76  mm  gun,  given  the  following  conditions:  -  1.654; 

s  -  0.4693;  /A-  18.44;  q  -  6.5;  pQ  -  30,000;  f  -  900,000;  a  -  1; 

$  -  1.6;  0  •  0.2;  w-  0.930;  -  1.05;  g  -  98.1,  using  strip  powder 

type  SP  (1  x  18  x  40  mm).  It  is  assumed  that  the  powder  burns  in 
conformity  with  the  physical  law  derived  from  the  pressure-bomb  test 
presented  in  Table  5  (correlation  between  I  andqO. 

To  start  with,  we  find: 

log  -  0.2970;  log  -  1.1709;  log  vjp  -  log  ^  )  -  8.0803; 


1  - 


-  2.284;  log  t  -  1.8710;  *  -  0.038;  log  —  -  6.2512. 

»  II  B 


F-TS-7327-RE 


659 


r 


O!-: 


i  si  J.  ICC. 


-Mttm 


*at 


From  the  value  found  we  compute  the  value  of  t^  by  inverse 
interpolation  with  the  aid  of  the  data  in  lable  2  (cf.  Example  2, 
Subsection  7) 

t0  -  0.020556  sec 

For  this  value  of  t^  we  determine  the  initial  impulse  I  by 
direct  Interpolation  with  the  aid  of  lable  4  (cf.  Example  1, 
Subsection  6) : 

Ig  -  76  kg.sec.d®**2 

Ye  find  the  velocity  of  the  projectile  at  the  end  of  burning 

s 

of  the  powder,  keeping  in  Bind  that  I  -  823  (cf.  Table  5): 

K 

v  -  ~(I„  -  1A)  -  5040  dm. sec”1 
K  ?m  X  0 

Bt  choosing  n  -  20  as  the  number  of  sections,  we  obtain  the  - 
step  as  being: 

VK  -1 

h  •  —  ■  252  dm* sec  1 

20 

so  that 

On  the  basis  of  Table  5,  we  plot  the  ty,  I)  curve  to  the 
following  scale:  for  4*,  1  mm  *  0.001;  for  I,  1  mm  *  1  kg* sec .dm"2. 
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We  read  on  it  the  vi.hics  of  the  gas  in!  low  v^*  corresponding  to  e.unl 


intervals 


823  -  76 
20 


37.35  kg*soc*dra_“ 


loi  the  pressure  impulse,  or,  what  is  the  same  thing,  to  ecjual 
intervals  h  -  252  dm»sec“*  lor  the  velocity  of  the  projectile. 

These  values  of  v^'are  given  subsequently  in  the  working  form 
used  for  performing  the  computations  (in  the  fifth  row). 
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It 


1 

V 

- 

0 

252 

0 

252 

5 

2 

2  log  v 

4.8028 

5. 4 

3 

-log  V“ 

Op 

9.91*7 

U . 9 

4 

j«4 

V0d _ 

1.122b 

• 

3.3 

5 

0.038 

0.092 

0.1 

2 

6 

V 

"I  *v 

/ 

-0.001 

-0.0 

1  v«p 

„2 

7 

*-  “f—  i 

0.091 

0.1 

V»P  / 

t 

log  fa 

■k, 

1.8710 

1.8710 

1.8 

9 

log  4* 

2.5793 

2.9638 

1 

1.2 

10 

log  Ia4* 

2.4503 

2.8348 

1.0 

11 

«A 

2.284 

2.284 
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The  working-  form  can  be  broken  down  into  three  main  sections. 

The  upper  section  contains  linos  ]-24  and  is  reserved  for  computations 
necessary  for  determining  the  logarithm  of  the  auxiliary  function 
In  accordance  with  the  relations  (*).  Consequently,  the  pressure 
Of  the  powder  gases  is  also  found  at  the  same  time.  The  values  for 
the  gas  inflow  in  line  5  are  first  read  from  the  ty,  I)  curve  in 
the  case  of  the  physical  law  of  burning  of -the  powder  or  computed 
in  accordance  with  the  corresponding  formula  in  the  case  of  the 
gsonetrlc  law  of  burning. 

The  upper  section  is  the  most  involved  part  of  the  work  and 
aakes  the  use  of  four-place  logarithms  obligatory. 

In  the  middle  section,  which  consists  of  lines  25-34,  are 
performed  computations  necessary  for  the  use  of  the  working  formula 
for  the  numerical  integration  of  a  first-order  equation: 


+  —  A2*.  «  +  —A3*  , 

12  i~2  s  i-3' 


i 

i 


and  for  finding  the  path  of  the  projectile:  j  +Afi_y 

Finally,  the  lower  section,  consisting  of  lines  35-47,  is  reserved 
for  finding  the  time  of  motion  of  the  projectile  in  the  bore  in  the 
form  of  a  definite  integral,  using  the  following  .working  formula: 


At,  -  F,  ♦  i-AF,  -  ±  A,  ♦ 


.  This  section  is  best  filled  after  completion  of  the  numerical 
integration  of  the  equation  for  the  path  of  the  projectile: 
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i.e.,  after  the  two  sections  above  it  are  filled. 

It  is  useful  to  point  out  some  of  the  characteristic  features 
of  the  computations. 

'  a)  The  path  of  the  projectile  is  determined  with  an  accuracy 
of  0.001  dm;  the  time  t,  within  0.00001  seconds. 

b)  To  facilitate  computations,  only  one  approximation  is 
made,  and  only  one  wedge  is  filled  out  (cf.  line  26,  first  column 
of  the  working  form). 

c)  At  first,  the  differences  of  the  auxiliary  function  up  to 
and  including  the  second  order  are  introduced. 

d)  The  third  differences  of  this  function  are  utilized  only 

after  p  is  passed. 

81 

e)  The  second  difference  of  the  function  $  appears  only  in 
the  second  column  of  the  normal  computations. 

f)  The  rows  which  are  not  filled  are:  lines  2,  3,  4,  6,  7,  21, 

22,  23,  24,  27,  28,  30,  31,  32,  and  35-47  in  the  first  column,  lines 
27-33  and  35-47  in  the  second  column  of  the  approximation,  lines  2-24, 

26,  27,  28,  31,  and  32  of  the  first  column  of  normal  computations, 
and,  finally,  lines  2-12,  27,  28,  31,  and  32  of  the  second  column  of 
normal  computations. 

g)  Line  34  is  the  leading  row  and  is  filled  first  in  each  column. 

For  a  more  definite  conception  of  the  character  of  the  work  Involved, 

when  filling  in  the  two  upper  sections  of  the  working  form  it  would 
be  desirable  to  describe  in  detail  the  computations  performed  for 
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OHO  of  trie  columns,  for  example,  for  the  column  of  normal  computations 
■^iul  *  — '  504. 

34)  The  ..umbers  ’n  lines  33  and  34  of  the  preceding  column  are 

added : 


»< 

l 


t  -  0.198  +  0.064  -  0.26?. 


since : 


lt  -  / 1-1  +  *lt_xs 


1)  Ve  add  the  integration  step  h: 

« t 

»  \ 

v  -  252  +  h  -  252  +  252  -  504; 

2)  The  number  in  the  preceding  row,  v  -  504,  gives: 

2  log  v  -  5.4048; 

3)  The  complement  of  the  logarithm  of  vjj^  is  taken: 

-log  v„  -  9.9197; 
nP 

4)  The  logarithms  in  lines  2  and  3  are  added: 

g  -  5.4048  +*9.9197  -*3.3245; 
v2 
np 

5)  The  value  of  ^  is  taken  from  the  {y,  I  )  curve t 

^  -  0.166; 
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6)  The  logarithm  in  line  4  is  used  to  find  the  number  vithin 
0.001: 

/ 

2 

—  -  0.002, 

vnP 

•  / 

7)  The  number  in  line  6  is  subtracted  from  the  number  in  line 

5: 

- -  0.166  -  0.002  -  0.164; 

V2 

Op 

8)  The  logarithm  of  l has  been  found  by  preliminary  computations: 

log  la  -  1.8710; 

9)  The  number  in  line  5  gives  the  logarithm: 

log  *  -  7.2201; 

10)  Ihe  logarithms  in  lines  8  and  9  above  are  added: 
log  la  1.8710  +7.2201  -  1.0911; 

11)  The  reduced  length  has  been  found  by  preliminary  computations: 

I  -  2.284; 

A 

12)  The  logarithm  in  line  10  is  used  to  find  the  number  vithin 

0.001: 
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la-t  "  0-123; 


13)  The  number  in  line  12  is  subtracted  from  the  number  in 


line  11: 


L  -L  -  l  *  «  2.284  -  0.123  -  2.161; 
T  A  a  1 


14j  The  number  in  the  leading  row,  line  34,  is  taken: 
/  I  -  0.262, 


15)  The  numbers  in  lines  13  and  14  are  added: 


l^+  t  -  2.161  +  0.262  -  2.423; 


16)  This  logarithm  has  been  obtained  by  preliminary  computations: 


log  —  -  6.2512; 

8 


17)  The  logarithm  of  the  number  in  line  7  is  found: 


N)- 


1.2148; 


18)  The  complement  of  the  logarithm  of  the  number  in  line  15 


is  taken: 


•log(iH  +  f)  -  -log  2.423  -  1.6257; 
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19)  The 'logarithms  in  lines  16,  17,  and  18  are  added: 


log  p  -  6.2512  +  1.2148  +  1.6257  -  5.0817; 


20)  The  logarithm  in  line  19  is  used  to  find  the  pressure 


at-  v  *  504 


p  «  1207  kg.cm" 


21)  This  logarithm  has  been  obtained  by  preliminary  computations: 


(”*)- 


1.5723; 


22)  Prom  the  number  in  line  1  we  have: 


log  v  -  2.7024 


23)  The  complement  of  the  logarithm  in  line  19  is  taken: 


-log  p  -  -5.0817  -  6.9183 


24)  The  logarithms  in  lines  21,  22,  and  23  are  added 


log  *  -  1.5723  +  2.7024  +  6.9183  -  1.1930. 


We  now  proceed  to  the  next,  middle,  section.  \ 

25)  The  logarithm  in  line  24  is  used  to  obtain  the  number  with 


an  accuracy  of  0.001: 


$2  -  0.156; 
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It  is  now  possible  to  take  the  next  step,  i.e.,  proceed  to 
the  next  column  with  v  -  756,  starting  the  computations  therein, 

.ms  always,  fro*  the  leading  row,  line  34. 

*  (  At  the  auzzle,/we  have:  -  583^m/sec;  p^  -  625  kg 'em2. 

13)  Numerical  Integration  of  the  Second-Order  Differential 

Equation. 

In  the  general  case,  a  differential  equation  of  the  second 
order  contains  the  following  components:  the  independent  variable  x, 
the  function  y  itself,  the  first  derivative  y'  of  the  function  y 
with  respect  to  x,  and  the  secopd  derivative  y"  of  the  sane  function 
with  respect  to  x.  Consequently,  this  equation  can  be  represented 
in  the  following  for*: 


t 


F(x,  y,  y*,  y")  -  0; 

where  the  synbol  T  represents  an  elementary  function  or  a  combination 
of  such  elementary  functions. 

Numerical  integration  *akes  it  possible  to  solve  this  equation 


in  any  for*,  provided  only  that  the  given  equation  permits  the 


deteralnatlon  of  the  second  derivative  as  an  explicit  function  of 


all  the  remaining  variable  quantities: 


y"  -  f(x,  y,  y*). 

/ 

♦  . 


In  this  numerical  Integration  the  increment  of  the  derivative 
Ay*  is  determined  by  the  aid  of  the  following  formula: 
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where  use  is  made  of  the  following  new  auxiliary  function: 


♦  -  hy”, 


expressed  in  terms  of  the  second  derivative  of  the  desired  function. 


It  is  not  necessary  lo  derive  the  fundamental  formula  for  Ay '  by 


means  of  y"  and  differences  of  this  second  derivative,  inasmuch 


as  this  formula  is  already  obtained  from  the  previously  derived 


relation 


Ayi  ■ h ( yl +  +  n*yi-2 4  IA-3 


by  simply  replacing  yi  by  y^  and  yf  by  y^’ 


On  the  other  hand  it  is  necessary  to  derive  a  working  formula 


for  the  increment  Ay  of  the  function  itself  in  terms  of  differences 


of  the  second  derivative  rather  than  of  the  first.  Such  a  formula 


should  simplify  the  work,  because  in  computing  Ay  it  will  be  possible 


to  make  use  of  the  already  available  differences  of  the  auxiliary 


function  and  it  will  not  be  necessary  to  find  the  differences 


of  still  another  auxiliary  function: 


F  -  hy* 
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For  this  purpose,  we  shall  utilize  the  equality 


-f 


AyA  -  l  y’dx. 


But  in  all  cases 


'  -  y ;  +  Ay*  -  y[  +  j* 


y"dx. 


Let  us  replace  here  the  derivative  y"  by  the  following  inter¬ 
polation  function: 

3  2 

y"  -  Aq(x  -  xA>  +  Aj(x  -  xi)  +  A2(x  -  x^  +  A3, 


and  then 


* 

'  "  yi  +  “  xi^3  +  Al<x  ”  XA>2  +  A2<x  ~  XA>  +  a3 _7  dx  - 


•  F*  +  A0(x  “  xi)4  +  *4  Al^x  ”  xi>3  +  \  A2^x  “  xi 


)2  +  A3(x  -  x j ) , 


so  that,  consequently 
*i+l  [ 

Ayi  "  j  I  yi  +  T  A0<x  “  Xi^4  *  T  Al^x  “  xi)3  + 


F-TS-7327-RE 


Of */*:•!  A.  MLY-  K  S>  L  J. 

W1WW.X ,  a**.** Sftv^rw  rynywpaa 


W*4 


. . n  ■»— wi 


1  2 

A9(x  -  x^)  4  A^(x  -  Xj) 


4  T  a2 


dx  - 


>i(x  -  xi>  +  jo  Aolx  '  V5  + 


l+i 


+  H  Vx  -  X1>4  +  T  Vx 


Xi)3  +  —  A  (x 
1  9  3 


Xi> 


whence 


A y  “  h 


y!  +  i-  Ah4  +  ~  A.h3  +  —  A_h2  +  — A.h 

i  20  0  12  1  6  2  2  3 


because 


xi+l  "  Xi  “  h* 


We  shall  substitute  into  the  last  expression  the  following  values 
for  the  coefficients: 

,2y”  „  +  A3y," 


A3  +1 

A  yi-3- 

0  61.3  ’ 


A1  - 


*  >i-2  +  a  yj'-3 

2h2 


A2  - 


CAy-C,!  4.  3A2y{’_2  +  2A3yj’_3 
6h 


I  ”  y j[% 


because  the  part  of  the  function  is  here  played  by  the  derivative 

1 

y",  which  has  been  replaced  by  the  interpolation  function. 

After  obvious  transformations  this  substitution  gives: 
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Ayt  -  h 


+  h 


*  e^i-i 


8 


i-2 


1  a3  .. 

—  A  V 

lu  yi-3 


if*  instead  of  the  iuconvenient  number  — — ■  we  take  the  closely 

180 

spproching  it  number: 


18  _  1 
180  10 


It  remains  necessary  to  make  use  once  again  of  the  auxiliary 

I  function: 

•  ( 

\ 

$  -  hy". 


and  then,  finally,  in  the  numerical  integration  of  the  second- 
order  differential  equation,  we  shall  have  to  deal  with  working 
formulas  of  the  following  type: 

Ayi  "  *i  +  TA*i-l  +  H  *2*i-2  +  \ A  3^*i-3  * 


Ay,  "  h 


(*i  +  +  7 


A  *1-1  +  7*1-2  +  ^a3*1-3, 


To  simplify  the  computations,  it  is  better,  however,  to  conduct 
the  numerical  integration  with  so  small  a  step  h  as  would  enable  us 
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to  do  without  third  differences. 

At  the  start  of  the  numerical  integration,  as  before,  most 
frequent  use  is  made  of  the  method  of  successive  approximations. 

The  number  of  rows  in  the  lower  part  of  the  form  will  increase 
(in  comparison  with  the  number  of  rows  in  the  lower  part  of  the  form 
lor  the  integration  of  the/f i rst-order  equation)  by  at  least  seven 
rows,  because  in  addition  to  the  ten  rows: 

♦  »  A  A  -~A*,  r|A~*i  —  A3$»  T,  Ay'  -  h£,  y', 

corresponding  to  the  lower  part  of  the  form  for  the  numerical 
integration  of  the  first-order  equation,  there  will  also  appear 
the  following  additional  rows: 

i*'  iA*-  i*2*-  ^A3*>ri’  Ay  ■  hEi’  y- 

In  the  method  under  consideration,  the  computations  themselves 
are  in  no  way  different  from  similar  computations  for  the  solution 
of  the  first-order  equation,  so  that  the  need  for  citing  a  special 
example  is  obviated. 

14)  Use  of  Numerical  Integration  of  Second-Order 
Equations  with  Argument  t. 

In  this  process  the  leading  part  is  played  by  the  equation  for 
the  forward  motion  of  the  projectile: 

dv 

tf>®  dt  "  ®P' 

from  which  it  follows  that 
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because 
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/ 


(107) 


« 


<  ■ 


if  the  phenoaenon  of  recoil  is  not  considered  in  the  explicit  form. 

Equation  (107)  is  subject  to  numerical  integration  by  means 
of  the  following  relations: 


P 


and 


•  -A 

which  have  already  been  'employed  earlier.  The  auxiliary  function  $ 
is  equal  in  this  case  to: 


* 


sh 

<pm 


where  the  step  h  is  already  a  certain  time  interval  of  the  motion 
of  the  projectile  in  the  bore.  This  step  is  chosen  in  advance. 
The  working  formulas  for  the  numerical  integration  will  be: 


1  - 6  o  3  3 

Av.  -  *4  +  — !  +  —  A  *4  o  +  — ‘  A 

1  t  2  12  a  i-3’ 


8 


A'i  -  li  ^  ♦  i.*t  ♦  A*l-1  +  Ta2*1-2  +  To*3*  1-3)' 


because  the  function  y  to  be  determined  represents  in  this  case  the 
relative  path  /  of  the  projectile  in  the  bore,  and  its  derivative  y' 
is  the  velocity  of  the  projectile  y. . 
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The  preliminary  computations  are  reduced  to  the  determination 
of  constants: 


1)  log  ■—;  2)  log  Sg;  3)  log  vjj^;  4)/^  ;  5)  log  (q;  6)  tVQ; 


/ 

/ 


7)  log 


as  in  the  case  of  numerical  integration  of  a  first-order  equation. 

The  values  of  I  ,  I *  and  v  are  found  in  an  analogous  manner. 

U  K  K. 

•v . 

The  gas  inflow  4>  is  read  off  the  (1,  v^»)  curve,  but,  in  contrast 
with  the  preceding  case,  it  is  necessary  here  to  find  for  each  point, 
during  the  process  of  integration  itself: 


1  -  I. 


ts, 

s 


because  the  integration  gives  values  of  the  velocity  v  which  are 
separated  from  one  another  by  unequal  intervals.  Once  the  value  of 
the  pressure  impulse  is  had,  we  can  obtain  the  required  value  of 
from  the  I,  curve.  Consequently,  this  curve  cannot  be  utilized 
in  advance,  but  it  must  be  available  in  the  course  of  the  entire  work 
of  numerical  integration. 

In  the  case  of  the  geometric  law  of  burning,  the  I,  curve  is 
replaced  by  the  following  relation: 


+  -  +0  + 


*e0?m  **<f2a2 


—  v  + 


si. 


s2!2 

8  R 


V2. 
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Fig.  153  -  Velocity,  Path,  and  Pressure  Curves. 

1)  Tine;  2)  pressure  curve;  3)  velocity  curve; 

4)  path  curve. 

Figure  153  contains  curves  for  p,  v,  and  l  as  functions  of  time 
t,  obtained  by  numerical  integration  with  respect  to  the  argument  t. 
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f  4) 
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log  * 
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6.2512 

6.251  t2f<si 

2 . 5793 

2.6532 

1.6466 

1.646-C 

4.4771 

4.5508 

3.2270 

3.2270 

2.258  2.271 


8.2512  6.2512  ' 

2.6532  2.7559 

1.6463  1.6438 


4.5509  4.6509 

3.2370  ‘3.2270 


■v 


Compu talon  Form  with  -Ar^ui^-nt  t. 


*»•  V>i-»  '.^V'^^SfEgf 


,0«  la 

log  ^ 


log  fa* 


A 


log  p 
log 


(s 


log  + 
p  kg*ci*2 


-V* 

12 

8 


0.038 


2.256 

J> _ 

2.256 

6.2512 


1.7041 

300 

51 


GsD 


0.045 

0.000 

0.045 


1.8710  1.8710 

2.5793  2.6532 

2.4503  2.5242 

2.284  2.284 

-0.028  -0.033 


2.251 

0.006 

— 

2.257 

6.25U 


2.5793  2.6532 

1.6466  1.646*.! 


4.5508 

3.2270 

1.7778 

355 

300 

"  60 

51 

9 

(2) 

(7) 

0.045  0.057 

0 . 000  0 .  uOo 

0.045  0.057 

1.8710  1 . 87 iO 

2. €532  2.7559 

2.5242  2.6269 

2.284  2.284 

-0.033  .0.042 

2.251  2.242  ' 

0.007  0.029 

2.258  2.271 

6.2512  6.2512 

2.6532  2.7559 

1.6463  1.6438 

4.5509  4.6509 

3.2370  3.2270 


1.7779 

355 

60 

9 

•  (7) 


1.8779 

448 

76 


300 

51 

(2) 

(7) 


21 


2.25f> 


2.257 


2.2  58 


2.27i 
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2.258  I  2.271 


6.2512  6.2512 


2.6532  2.7559 


1.6463  1.6438 


2.294 


i.  .  1  DO 


0.435  \ 


2.330  2.395 


6.2512  6. -.>12  6.2512 


2.8976  1.0719  *1.2304 


1.6394 


l.*.o26  1.6209 


4.7882  4.9557  5.1025 

3.2270  3.2270  3.2270 


2.1827 


0.007  0.0022 


0.000  0.007 


0.007  0.022  0.040  0.066 


0.000  0.007  0.029  0.069 


2.3295 


0.159 


0.135  I  0.239 


CHAPTER  2.  SOLU'l  ION  BY  EXPANSION  IN  TAYLOR* S  SERI KS 

V 

Integration  of  equations  of  internal  ballistics  describing  the 
relations  existing  between^  the  fundamental  elements  of  a  shot  leads 

/ 

to  rather  complex  integrals,  which  can  be  solved  only  by  means  of  y 
quadratures  with  any  desired  degree  of  accuracy  (Professor  Drozdov's 
solution).  In  this  connection,  in  order  to  make  integration  possible, 
some  of  the  variable  parameters  in  the  fundamental  equations  (0,  u^, 

Pq>  etc.)  are  usually  assumed  to  be  constant. 

The  method  of  numerical  integration  makes  it  possible  to  solve 

a  system  of  differential  equations  not  only  without  simplifying  the 

« 

functions  entering  same,  but  even  by  assigning  variable  values  to 

/ 

those  quantities  which  are  usually  assumed  to  be  constant.  This 
makes  it  possible  to  solve  the  problem  on  the  basis  of  any  desired 
hypotheses  concerning  the  character  of  burning  of  the  powder,  the 
law  of  resistance  to  motion,  the  design  of  the  bore,  etc. 

In  solving  the  problem  by  the  method  of  numerical  integration, 
it  is  necessary  to  proceed  successively  from  one  value  of  the 
argument  to  another  by  the  addition  of  its  finite  differences,  starting 
at  the  very  beginning.  For  this  reason,  for  example,  it  is  not 
possible  to  determine  in  advance  the  values  of  pm  or  the  values  of 
the  variables  v,  (  and  p  corresponding  to  the  end  of  burning  of  the 
powder,  it  being  nec <  ary,  instead,  to  compute  successively,  point 
by  point,  the  elements  of  the  curves  of  pressure  p,  the  path  of  the 
projectile  I,  its  velocity  v,  etc.  This  constitutes  the  disadvantage 
of  this  method.  Moreover,  the  method  of  numerical  integration  gives 
the  relation  between  the  individual  variables  only  in  the  form  of 
numerical  tables,  rather  than  in  the  form  of  analytical  formulas. 
F-TS-7327-RE  680 


^  Mrvrn  *>  «•  r\ 

. 


In  spite  of  these  disadvantages,  the*  method  of  numerical 
integration  rn.ty  serve  .is  a  means  for  indirect  verification  of  the 
degree  of  accuracy  obtained  by  the  aid  of  the  various  approximate 
analytical  methods  available. 

4 

In  this  connection,  when  developing  new  theoretical  piobleos, 
numerical  integration  may  be  utilized  for  determining  the  errors 
involved  in  the  transition  from  exact  equations  and  relations  to 
others  that  are  less  exact  but  more  convenient  from  the  analytical 
point  of  view.  Numerical  integration  may  be  employed  with  equal 
success  both  in  the  case  of  the  geometric  law  of  burning  and  in 
the  case  of  the  more  complex  physical  law  of  burning;  it  may  also 
be  applied  to  barrels  having  a  bore  of  variable  cross  section. 

In  the  USSR  the  method  of  numerical  integration  was  first 
applied  to  the  solution  of  ballistic  problems  by  V.M.  Trofimov  in 
1918.  This  method  has  been  developed  in  great  detail  by  G.V. 

Oppokov,  who  employed  the  method  of  finite  differences  (1931-1938) 
discussed  above. 

» 

In  1932  P.V.  Melentyev  proposed  to  apply  the  method  of  expansion 
in  laylor's  series  for  the  numerical  solution  of  equations  in 
ballistics,  and  this  method,  after  being  subjected  to  certain 
modifications,  has  been  found  to  be  sufficiently  convenient. 

Investigation  of  the  fundamental  relations  expressing  the 
conditions  accompanying  a  shot  shows  that  all  the  principal  elements 
(/,  v,  *V,  and  p)  can  be  expressed  in  one  form  or  another  as  functions 
of  the  path  I  and  of  its  derivatives  with  respect  to  time  up  to  the 
third  order  inclusive. 


As  a  natter  of  fact,  taking  the  fundamental  system  of  equations 
expressing  the  relationship  between  the  elements  of  a  shot,  we  obtain 
the  following. 

1)  The  fundamental  equation  of  pyrodynamics  or  the  equation 
of  transformation  of  energy: 


PS<f  +  I)  “  f«*  -  —  ‘fmv2. 


2)  The  equations  expressing  the  law  of  burning  of  the  powder: 


-  Kz(l  +  Xz)  ; 


de  m 


u  -  u,p; 


dq»  1  k 

-i-  -  - 6U.  P  -  — —  6u  D. 

dt  Ai  1  el  1 


3)  Ihe  equation  of  motion: 


ps  -  «fm  -  vfm  m  I  ” , 


where 


dt  t. 

w  -  —  -  r, 

dt  t 


v  being  related  to  z  by  the  following  equation: 


v  -  -  (z  -  z  ), 

«f»  0 
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whence 


«f  m 

7  “  zo  "  —  V< 

si* 

K 


All  the  variables  entering  into  the  fundamental  system  of  equations 
can  be  expressed  in  terms  of  the  path  /  and  of  its  derivatives,  since 


v  *  K>  7  "  z0  + 


I’,  -  being  a  function  of  z  -  will  also  be 


expressed  as  a  function  of  t^,  the  pressure  is  proportional  to  f*^, 

dp  , 

and  the  derivative  ~  is  proportional  to  l .  Consequently,  if  the 

time  t  of  travel  of  the  projectile  through  the  bore  is  taken  as  the 

independent  variable,  and  the  path  of  the  projectile  l  is  taken  as 

the  function  to  be  expanded  in  a  series,  it  becomes  possible  to 

employ  Taylor’s  series  for  finding  the  value  of  the  path  and 

of  its  derivatives  for  the  neighboring  segment  corresponding  to  the 

time  tn+i  *  t  +  At  -  tn  +  h,  provided  that  the  values  of  the  path 

/  and  of  its  derivatives  for  the  preceding  instant  t  are  known, 
n  n 

It  is  thus  possible  to  find  all  the  elements  of  burning  of  the  powder 
and  of  the  motion  of  the  projectile  during  a  shot,  i.e.,  z,  v,  p, 

( ,  and  t. 

Let  it  be  assumed  that  for  a  certain  instant  of  time  tn  the  path 

/„  and  its  derivatives  with  respect  to  time  are  known; 

if  a  sufficiently  small  increment  of  time  At  ■  h  is  assumed  and 

consideration  is  limited  to  derivatives  up  to  and  including  the  third 

order,  then,  in  accordance  with  Taylor's  formula,  we  shall  have  for 

t  -  t_  +  h: 
n+1  n 
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_  -  5  ******  2* 


n4l 


-I  4  h/ 


^  / 


n 


4  HI  /"’ 

n  2*3  n 


(108) 


i' 

i 


Differentiating  with  respect  to  t  and  rejecting  the  terms 
1 V 

containing  (*T,  i.e.,  considering  that  t is  constant  over  the 
given  interval  At  and  equals  its  mean  value*  we  obtain: 


I*  ,  -  t '  4  hi" 
•n4l  n  n 


2 

4  /,M  ; 

/  n 


/"  -  I”  +  hi"’  “  lB  +  h 

*n41  n  n  av.  n 


I  •••  +  f  •*  * 

‘n  +ln4l 


n  av. 


where 


r’4 


2 


®-iis  the  mean  value  of  the  third  derivative  in  the  interval 


under  consideration  (Fig.  154).  From  the  last  equation  we  obtain  the 
following  value  for  t"',:  ^n4l  •  ~  *n4-l  -  ~  l'n  -  in  . 


Fig.  154  -  Diagram  for  ("  and  J”* . 

As  has  Men  shown  by  P.V.  Melentyev,  it  is  more  convenient  to 

compute  not  the  derivatives  themselves*  but*  rather*  the  quantities 

2  h^ 

proportional  to  them,  namely:  hi',  hi"  and  —  Therefore*  by 

i.3 

multiplying  by  h  and  by  2- *  we  will  obtain  the  following 


equations: 
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'  h|r  -  h|'.+  h2!'*  *  ~  1  ; 

n«l  it  n  2  n 


(109) 


12 _  i 

2  h+1 


h2  j " 
ln+l 


h2  | " 


r  i  mi 
T  * 

*  n 


(110) 


Comparing  these  two  equations  with  the  initial  equation  (108), 

it  will  be  seen  that  f  enters  everywhere  without  a  coefficient,  /* 

enters  with  the  coefficient  h,  /”  enters  with  the  coefficient  h2, 

h3 

and  /"*  enters  with  the  coefficient  This  considerably  accelerates 

the  subsequent,  computations. 

/ 

The  expression  for  the  second  derivative  / "  will  be: 

r:  -  —  p 


or,  multiplying  by  hz : 


h2l"  -  h2  £-p. 
t  q»m 


(111) 


L,  combining  the  resulting  values  for  the  path  f  and  its 
derivatives  with  the  equations  of  the  fundamental  system,  we  obtain 
the  totality  of  formulas  necessary  for  the  solution  in  the  following 
form  and  sequence,  which  corresponds  to  the  order  of  their  application, 
the  constants  encountered  being  designated  below  as  follows: 


Jts  -  k., 

SiK 


■  *2’  IS  ■  V  T*  -  V  T  -  t) 


■v  a. 
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hr  »  hl'  4  h2r  ♦  —  I'"  ; 

n + 1  n  ~ 


n+1 


n  2  n 


hi*  , 
n+1 


System  I 


z»+i  ■  7o +  rr  \*i m  'u  +  kiv,.+i: 


+.+1  m*\+l  +  *Xzn+l: 


3 

1  ,  -  I  +  hi'  +  4  »»2I"  +  —  ; 

n+1  n  n  2  n  i  ?  '«  ’ 


^n+1 


v2 

vn+l 


r«P 


n+1 


8  I  .  / 

wn+1  ‘n+1 


n  32 


.  k  Vl  ~  V-+1 

2  <4-  »+n+l  +  <„+l 


2  |  n  _  ».2  S  _  t2. 


-2'n"+i  '  » 


<f  a 


Pn+1  '  h  V„+l ; 


1 

£  , 
a  *  - 


2  ’  n+1 


■  h2/»*  -  h2 f  "  -  tL,  l,u 

*  n+1  ft  ‘n  2  'n  * 


1 


(1) 


(II) 


(III) 


(IV) 


(V) 


(VI) 


(VII) 


(VIII) 


The  subscript  (n+1)  designates  those  values  of  the  derivatives 
at  the  end  of  the  given  interval  of  tine  which,  in  the  process  of 
computation,  are  transferred  from  the  column  being  computed  into  the 
corresponding  rows  of  the  right-hand  neighboring  column;  however, 
the  transferred  values  now  bear  the  index  n  because  they  characterize 
the  initial  value  of  the  given  quantity  in  the  next  column. 
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In  order  to  perform  t fie  computation  by  means  of  this  totality 
of  formulas,  the  values  of  /  and  of  its  derivatives  at  the  start 
of  the  projectile's  motion,  i.e.,  at  the  instant  t  -  0,  must  bo 
known.  Since  at  the  start  of  motion  the  path  (  and  the  speed  v 
are  equal  to  zero,  we  obtain: 

(I) 0  -  0,  (|’)0  -  (v)Q  -  0,  ,<l">0  -  k4p0,  h2(|")0  -  h2k4P(), 


where  Pq  is  the  forcing  pressure  usually  specified  beforehand.  As 
regards  the  third  derivative  (1'")^,  we  shall  first  find  an  expression 
for  it  at  the  present  instant  in  the  form  of  . 

To  determine  ,  we  differentiate  the  equation  /*'  -  k4p  with 


respect  to  t: 


*;■  -  v; 


But  the  quantity  p^  has  already  been  derived: 


-  v ( 1  +  0)  , 


for  the  start  of  motion  when  1-0,  v  -  0,  p  -  p.  ,  and  ;  we 


will  therefore  obtain: 


*  \5s_ 

"I  4, 


"  v  I  9  r  M  *  '()'  i 


quantity  1  +  77-  -  7=-,  and  this  expression  is  therefore 

1 
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sometimes  given  in  the  following  form: 


X6„  / 

0  'A 


(pt}0  "  k2  fjTpU* 

K  *0 


Consequently,  the  value  of  the  third  derivative  for  the  initial 
Instant  is  likewise  known: 


,so  '4 


£_/•««..  h_v  (n<)  .  L  v  k  2 

2*0  2  Wo  2  V2  1K 


K 


and  it  is  possible  to  begin  the  successive  solution  of  System  (I) 
first  for  the  first  column  corresponding  to  the  interval  of  time 
At  *  h,  then  for  the  second  column,  etc.,  thus  obtaining  a  successive 
series  of  values  for  /,  v,  z,  ^  and  p  as  functions  of  t. 

The  quantity  h  *  At  must  be  so  chosen  as  to  obtain  10-15  columns 
for  the  period  of  burning  of  the  powder,  which  will  give  a 
corresponding  number  of  points  for  each  of  the  quantities  p,  v,  /, 
and 

Since  the  time  of  burning  is  fundamentally  determined  by  the 
thickness  of  the  powder,  the  interval  of  time  h  -  At  may  be  taken 
approximately  according  to  the  formula: 


h  a*  0.001  ej. 


where  e  ■■  —  the  thickness  of  the  powder  in  millimeters,  h  being 
1  * 

rounded  off  to  one  or  two  significant  figures  (to  5  in  the  second 
significant  figure).  For  example,  if  2e^  -  1.28  mm: 
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h  -  O.Oul  (U.i>4)  -  O.UIHI64  -  O.OUOti  or  0.0i»l‘65 

Si  net*  vA  and  /  lire  known  n  t  least  approx  imn  le  I  v  in  advance, 

it  is  possible,  after  computing  tin*  average  lime  ot  motion  of  the 

2  / 

nroiectilo  t  *  - ,  to  take  lor  the;  value  of  the  time  step  (increment) 

'  ‘  t  iiv.  vA 

At  -  h  ~  ,  rounded  off  to  two  significant  figures. 

1  o 

Sequence  of  Computation.  All  the  constants  are  computed  first: 

1  ~  V  1  1 


X,  X,  X  X,  A,  4’y 


JL  +  a  -  1. 


f  +  4- 


©0  “\  +  4 


0 


Po 


2*, 


z  - 

0 


0 


*(<50  +  1) 


JK  “ 


_ 1 .  /  - 


s 


■  7  fa  -  ^  \  -'a  -V  f  -  5571*  h  *  1  «e 


1  -  —I  k2  ■  Sr;  k3  ’  IS  ■  ~r  <SB:,li  q»A»tity> ;  k_,  -  -±;S 

sl„ 


np 


h2,o 


2  h3  „„  h3  «*0  'a 

hVo=  — ■  —  ktk2  tk  —  po 


^0 
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The  sequence  oi  compulation  is  not  ntfectcd  regardless  ot  whether 

% 

•  •' 

the  computation  is  performed  for  a  degressive  or  a  progressive 
/  4  f 

powder. 

In  computing  the  segment  after  theytlecomj  osi  t  ion  of  the 
progressive  powder, it  is  necessary  to  substitute  for  the  usual 
formula  the  following  previously  derived  formula: 


'V-  % 


+  K  (z 

z 


l)/~  1  +  \2u 


-D7- 


x2(z  -  l) 


where  z  varies  from  1  to  1  +  ~ ,  and  rt0  and  X0  are  characteristics 
of  the  powder  form  after  decomposition. 

A  form  for  conducting  such  computations  is  presented  on  pages 
691-692. 


Table  7-c  -  Computation  Form  for  the  Solution  of  Problems  in 
ballistics  by  l.xjansion  in  'laylor's  Scries. 


Subse r i j  t  n 


Compu t  a  t  ion  Formu  1  .■  s 


h  -  0.0008 


hr  -  b i • 
n  +  l  n 


9  h 

h“|"  4 

n  2f« 


hf"  -  0 
0 


hi"  -  0.0958 
0 


—  -  0.0364 

2  0 


1 

2 

3 

0.0008 

0.0016 

0.0024 

0 

1.1322 

0.3856 

0.0958 

0.1928 

0.3546 

0.0364 

0.0606 

i 

0. 1012 

0.1322 

0.3856 

0.8414 

165.3 

482 

0.0001864 


klvn+l 


0.0308 

0.0297 


1  +  *^Zn  1 

1  n»  1 


-  0.087030 


2  -  2,850,000 


7  n  >  1 

0.0605 

*'.,i 

0.0641 

4- 

0.0002 

0.0639 

k3v^l 

0.0002 

^n+1  -  k3vLl 

0.0637 

k2(t„,l  -  k3v?al)  ‘ 

181500 

“  An4l 

r -T  S-7327-RK 


7-t  (Cunt fd. ; 


u  t  a  t  ioi  Y o rn.ij  1  a  s 


jVl  ~  V  1,-1 

?  «+„.!  "  '..I 


0.0008  0.0016  0.0024 


0  0.0600  0.3088 

0  0.1322  0.3856 


Li  -  +  + 

+  -jt>2r  +  -i-h3/”* 

^  n  3  n 


0.0479  0.0964  0.1773 


(l)0  -  0:  w>  -  0 


0.0121  0.0202  0.0337 


-T(--i)- 


1.065  21  /n+1  +  (A 


22  »+n<l 


0.0531 192 


23  »  ,  -  /  .  +  / 

n+1  n+1 


-  aVl 

- — r 


n+1,  ? 

Pn+1  “  “k«  cm 

B„ .  1 


h2 1 "  . 
n  +  1 


-h2|» 

n 


.  b!  /». 

2  n 


0.0600 

0.3088 

3.016 

3.016 

3.076 

3.325 

v  J68 

i 

0.  134 

3.008 

3.  191 

604 

1100 

0.  1928 

0.3546 

-0.0958 

-0. 1928 

-0.0364 

-0.0606 

2  ‘n+! 


0.0606 


0. 1012 


0  0.1322  0.3856 


0.0479  0.0964  0.1773 


4-  — I’"  0.0121  0.0202  0.0337 

3  2° 


0.0600  0.3088  0.9054 


*n+l  + 


•‘♦'n+l 


23  B  .  *  /  .  / 

n+1  u+1  A 


-  *'Vr 


3.016  3.016 


3.076  3.325 


i.:/J68  0.134 


3.008  3.191 


1100 


3.016 


From  lino  19  of 
preceding  column 


From  line  4  of  preceding 
column. 


—  (of  line  2  of  the 
2 

given  column) 


—  (of  line  3  of  the  given 
column) 


Into  line  15  of  next 
column 

In  all  columns 


-h 2 1" 
n 


0.1928  0.3546 


0.0958  -0.1928  -0.3546 


Into  lines  2  and  26  of 
next  column 

From  line  25  of  preceding 
column 


.  bf  i- 

2  n 


1-0.0364  -0.0606  -0.1012 


From  line  28  of  preceding 
column 


h^  ...» 

2  *n+l 


0.0606  0.1012 


Into  lines  3  and  27  of 
next  column 


I 


—  (of  lii.e  3  ol  the  eiwii 
3 


column) 


0.9054 


Into  line  15  of  next 
column 


i  * 

3.016  3.016  3.016  |  In  all  columns 


*  *n+l  +  *A 


3.076  3.325 


i'J6  £  0.134 


Bn+1  "  *n+l  +  (ft 


3.008  3.191 


’  -  *Vl 


604 


1100 


h 2(”  . 
n+1 


0.1928  0.3546 


•b2/' 


0.0958  -0.1928  -0.3546 


\  V  -  h-  /'” 
\  2  n 


0.0364  -0.0606  -0.1012 


t  h? 


M  ftfl 

2  'n+1 


0.0606  0.1012 


1  Formulas  for  the  second  period: 


f.fiiiiiV « 

Mr ) 


“  k3v2) ;  p  -  pK 


(0)' 


when  /-I, 


v  -  v„ 


P  -  P, 


are  determined  from  firsr  period  when  ^  -  1, 


Into  lines  2  and  26  of 
next  column 


From  line  25  of  preceding 
column 


From  line  28  of  precedi:. g 
column 


Into  lineF  3  and  27  of 
next  column 


yi'-i?- 


CKA.’.Ll '  A.  '*  r«  *  *• 


1 


The  extreme  left  column  contains  the  "computation  formulas" 
and  constants  of  System  (I);  in  the  next  column  to  the  right  these 
lormulas  are  broken  down  into  individual  operations,  which- are 
followed  in  the  computations. 

To  start  with,  the  first  column  (No.  1)  corresponding  to  the 
time  interval  0  to  h  is  filled  in  first.  In  this  column  the^ubscript 
n  relates  to  the  start  of  the  interval,  and  the  subscript  n+1  relates 
to  its  end;  for  this  column  n  -  0  and  n+1  ■  1. 

For  the  next  (second)  column,  n  -  1,  n+1  -  2,  etc.  For  the 

first  column,  computation  of  the  constants  gives  us  at  n  -  0  l  -  0 

n 

(the  path  at  the  start  of  the  motion),  which  we  write  on  line  15; 

h I’  •  0  (the  velocity  at  the  start  of  the  motion)  is  written  on 
n  2 

lines  1  and  16;  h  t"  -  h^l"  •  0.0958  is  written  on  lines  2  and  26; 

3  n  0 

—  is  written  on  lines  3  and  27.  Line  17  is  filled  with  — (h^l"), 

2  0  3  2  0 

and  line  18  with  )•  Thus,  all  the  quantities  with  the  subscript 

o  2  0 

n  •  0  are  inserted  in  the  first  column.  We  now  subject  them, to  the 

necessary  operations.  The  sum  of  the  first  th'ree  rows  gives  the 

fourth  h I’  ,  -  hi'  .  which  is  immediately  transferred  to  lines  1 
n+l  0+1 

and  16  of  the  neighboring  column  wherein,  provided  with  the  subscript 

n,  it  characterizes  the  value  of  this  quantity  at  the  start  of  the 

next  interval;  we  then  determine  vn+^,  7n+j»  H'n+l'  and  An+1* 

adding  the  four  rows  from  line  15  through  line  18,  we  obtain  in  line 

19  /  ,  -  the  path  of  the  projectile  at  the  end  of  the  given  Interval 

11  +  1 

of  time  -  and  this  quantity,  provided  with  the  subscript  n,  is 

transferred  to  line  15  of  the  neighboring  column.  After  determining 

pn+i  ■  and  multiplying  it  by  k^h^,  we  obtain  h2/jj^,  which  we 

®n+l 

write  on  line  25  of  the  first  column  and  on  lines  26  and  2  of  the 
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next  column,  where  this  quantity  acquires  the  subscript  n,  as  does  < 

1  2  ' 
also  -“(h  I"  ),  which  is  written. in  line  17  of  the  next  column. 

2  n+1 

After  performing  the  operations  indicated  in  the  form  with  lines 

3 

25,  26,  and  27,  we  obtain  in  line  28  of  the  first  column  k-/"’  , 

2  n-i  1 

which  we  immediately  transfer  to  lines  3  and  27  of  the  next  column, 

v  while  )  is  written  in  line  18  of  the  same  column. 

3  2  n+1 

Thus,  all  operations  with  the  quantities-  bearing  the  subscript 

n  in  the  second  column  arc  already  prepared,  and  the  second  column 

is  then  treated  In  the  same  manner  as  was  the  first. 

Constants  such  as  z_  in  line  7  and  /  in  line  20  are  Inserted 

U  A 

in  the  series  of  columns  in  advance. 

By  applying  the  same  rules  to  the  neighboring  second  column, 
we  shall  gradually,  step  by  step,  obtain  values  for  v,  z,  t, 
and  p  as  functions  of  t  -  (n+l)h,  and  this  is  continued  to  the  end 
f  of  burning  or  to  the  instant  of  emergence  of  the  projectile  from 

the  bore,  it  being  necessary  to  use  1  after  the  end  of  burning. 

In  performing  the  computations  it  is  necessary  to  exercise 
extreme  care  not  to  commit  any  errors,  because  an  error  in  one  of 
the  preceding  columns  will  distort  the  results  obtained  in  the 
succeeding  columns. 

It  is  best  to  follow  up  the  computation  of  the  data  in  each 
column  by  plotting  them  on  graph  paper  as  a  function  of  t..  In  so 
doing  an  error  in  the  given  column  will  cause  a  deviation  from  the 
regular  disposition  of  the  points  derived  from  the  preceding  columns, 

{  and  such  an  error  can  be  detected  and  corrected. 

As  a  criterion  of  accuracy,  it  is  also  useful  to  plot  on  the 

h3 

diagram  the  third  derivative  (or  9-  f'"  in  the  last  row),  which  *-■ 

2  n+1 
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should  first  increase,  thru  pass  through  a  maximum,  then  become 
zero  (p’  «  0)  at  the  instant  p|n  is  attained, and  thereupon  acquire- 
a  negative  value,  f Juctuat ing  s 1 ight 1 y  in  either  direction, 

h3 

The  instant  of  tim*'  cut  off  on  the  diagram  at  p*  -  0  or  — l"'  -  0 

/  1  2 

corresponds  to  the  irfstant  of  maximum  pressure,  and  all  elements 
for  it  are  best  taken  from  the  diagram 

The  time  t^  corresponding  to  the  end  of  burning  of  the  powder 

t 

is  determined  from  the  diagram  on  the  basis  of  the  vy,  t  curve  at 

l 

4*  •  1;  thereupon,  the  elements  corresponding  to  the  end  of  burning 
of  the  powder  for  this  time  are  found  by  interpolation.  If,  without 
changing  the  segments  At  -  h,  the  second  period  is  computed  as  a 
direct  continuation  of  the  first,  assuming  -  1  and  throughout, 

the  third  derivative  /"*  usually  begins  to  fluctuate,  sometimes 
entering  the  region  of  positive  values,  which  contradicts  the 
physical  nature  of  the  process  of  pressure  change. 

For  this  reason,  once  the  elements  of  the  end  of  burning 
VK’  ant*  pg  have  been  obtained  from  the  computation  of  the 
first  period,  tho  procedure  is  continued  by  adopting  the  same  step 
At  -  h  with  t£  as  the  starting  point  by  first  computing  the  values 
of  the  path  and  of  its  derivatives  for  the  start  of  the  second  period 
in  accordance  with  the  following  formulas: 


/(0)  -  fK;  hl’rni  "  hv*;  "  h2~P^  -  hVp*; 


<fm*K 


3  h3  (1  +  «)vKp 

h  ini  B  *} _  _S_  k 

2  (0)  ”2  <fm  (  +  I  ) 


—  -  -  —  k  (1  +  e )  ; 

*>  24  (>+,x 
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* 

V' 


f 


\t 


whereupon  they  are  written  in  the  initial  column  for  compute  ting 
the  data  of  the  second # period  and  subjected  to  the  sane  operations 
as  in  the  first  period,  with  the  sole  exceptions  that  ^  1  is 
assumed  in  line  11  and  that  lines  6-10  are  omitted. 

The  computation  is  continued  in  this  manner  until  I  ,  ^ 

If  t  »  /  ,  the  remaining  elements  v_  and  p_  are  obtained 

n+1  A  rS 

automatically  in  the  same  column  for  the  subscript  n+1 ;  if  l  n+1  >L’ 

A 

the  computations  in  this  column  must  be  carried  as  far  as  line  24, 
with  lines  6-10  omitted,  whereupon  the  value  of  /A  is  used  to  obtain 
tA  by  interpolation  in  the  last  column  for  the  purpose  of  subsequently 
obtaining  the  elements  pA  and  vA. 

Instead  of  expanding  in  a  series  after  obtaining  the  elements 
corresponding  to  the  end  of  burning  of  the  powder,  it  is  possible  to 
compute  v.  and  p  by  means  of  the  usual  second-period  formulas,  but 
without  determining  the  time  t. 

The  solution  by  expansion  in  a  series  is  applicable  to  both 
the  geometric  and  the  physical  law  of  burning  of  powder.  In  the 
latter  case: 


•■dt  -  V* 


from  which  we  have  the  following  expression  for  1: 


1  -  >0  4  ’S  *• 


and  the  correlation  between  ^  and  z  is  replaced  by  the  graphical 
dependence  of  I  upon  v^,  which  is  found  from  the  bomb  test. 
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The*  equation  for  is  replaced  by  the?  following  expression; 


<^>0  “  >:2^0  1 2  p0 

*0 


and : 


.3  h3  L 

£2 _ I  *"  ■  L  k  k  T  n 

2o  2  ^  2  0  ,2^  0> 

where  is  the  value  of  the  experimental  function  T  corresponding 
to  the  quantity  ^ 

In  the  case  of  ballooning  powders,  is  greater  than  in  the 
case  of  the  geometric  law  of  burning,  and  therefore  the  values  of 
both  the  first  derivative  p’  and  of  p  itself  will  increase  more 
rapidly,  and  the  maximum  pressure  will  occur  earlier.  If  the 
propellant  force  of  the  powder  is  the  same,  the  maximum  pressure  will 
be  greater  in  the  case  of  the  physical  law  of  burning  with  ballooning 
than  in  the  case  of  the  geometric  law  of  burning. 
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